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RECURSIVE COMPUTATION OF THE INVARIANT MEASURE OF 
A STOCHASTIC DIFFERENTIAL EQUATION DRIVEN BY 
A LEVY PROCESS 

By Fabien Panloup 

Universite Paris 6 

We study some recursive procedures based on exact or approx- 
imate Euler schemes with decreasing step to compute the invariant 
measure of Levy driven SDEs. We prove the convergence of these 
procedures toward the invariant measure under weak conditions on 
the moment of the Levy process and on the mean-reverting of the dy- 
namical system. We also show that an a.s. CLT for stable processes 
can be derived from our main results. Finally, we illustrate our results 
by several simulations. 

1. Introduction. 

1.1. Objectives and motivations. This paper is devoted to the computa- 
tion of the invariant measure (denoted by u) of ergodic stochastic processes 
which obey a stochastic differential equation (SDE) driven by a Levy pro- 
cess. Practically, we want to construct a sequence of empirical measures 
(p n (uj, dx)) n >i which can be recursively simulated and such that v n {uj,f) — > 
v{f) a.s. for a range of functions / containing bounded continuous functions. 

In the case of Brownian diffusions, some methods have already been de- 
veloped by several authors to approximate the invariant measure (see Sec- 
tion 1.3), but this paper seems to be the first one that deals with this problem 
in the case of general Levy driven SDEs. The motivation for this general- 
ization is the study of dynamical systems that are widely used in modeling. 
Indeed, there are many situations where the noise of the dynamical system 
is discontinuous or too intensive to be modeled by a Brownian motion. Let 
us consider an example that comes from the fragmentation-coalescence the- 
ory. In situations such as polymerization phenomenons, when temperature 
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is near to its critical value, molecules constantly break-up and recombine. 
This situation has been modeled by Berestycki [2] through what he terms 
EFC (Exchangeable Fragmentation-Coalescence) process. The mass of the 
dust generated by this process (see [2] for more details) is a solution to a 
mean-reverting SDE for which the noise component is driven by a subordi- 
nator (an increasing Levy process). We come back to this example in Section 
7. 

For other examples of situations where models that use a Levy driven 
SDE are adapted, we refer to Barndorff-Nielsen et al. [1] for examples in 
financial modeling (where ergodic Levy driven SDEs are usually used to 
model the volatility of a financial market), Protter-Talay [23] for examples in 
finance, electrical engineering, ... or Deng [8], who models the spot prices of 
electricity by a mean-reverting Brownian diffusion perturbed by a compound 
Poisson noise. 

1.2. The stochastic differential equation. According to the Levy-Khintchine 
decomposition (for this result and for introduction to Levy processes, see, 
e.g., Bertoin [4], Protter [22] or Sato [25]), an Revalued Levy process (L t ) 
with Levy measure tt admits the following decomposition: Lt = at + y/Q Wt + 
Yt + Nt, where a G M. 1 , Q is a symmetric positive / x I real matrix, (Wt) is 
a /-dimensional standard Brownian motion, (Yt) is a centered /-dimensional 
Levy process with jumps bounded by 1 and characteristic function given for 
every t > by 

t( f e^ -l-i(u,y)7T(dy)) 

and (Nt) is a compound Poisson process with parameters A = ir(\y\ > 1) 
and n(dy) = l/|yi > n7r(dy)/7r(|y| > 1) (A denotes the parameter for the wait- 
ing time between the jumps of ./V and n, the distribution of the jumps). 
Moreover, (Wt), (Y t ) and (Nt) are independent Levy processes. 

Following this decomposition, we consider an Revalued cdlg process (Xt) 
solution to the SDE 

(1) dX t = b(X t - ) dt + a(X t - ) dW t + n(X t - ) dZ t , 

where 6:R d i-^R d , a : M. d i-> M^ j (set of d x / real matrices) and n:R d i-^ 
M^z are continuous with sublinear growth and (Z%) is the sum of the jump 
components of the Levy process: Zt = Yt + Nf. 

In most papers dealing with Levy driven SDEs, the SDE reads dXt = 
f(X t ~) dL t , where (L t )t>o is a Levy process. Here, we separate each part of 
the Levy process because they act differently on the dynamical system. We 
isolate the drift term because it usually produces the mean-reverting effect 
(which in turn induces the ergodicity of the SDE) . The two other terms are 
both noises, but we distinguish them because they do not have the same 
behavior. 



E{e i{u ' Yt) } = exp 
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Remark 1. In (1) we chose to write the jump component by compen- 
sating the jumps smaller than 1, but it is obvious that, for every h > 0, {Xt) 
is also solution to 

(2) dX t = b h {X t - ) dt + a(X t - ) dW t + n(X t - ) dZ? 

with b h = b + J { \ y \ e(lM} VTr(dy) if h > 1, b h = b - J { \ y \ e{hjl]} y^{dy) if h < 1, 
and Zt = Y t h + , where the characteristic function of Y t h is given for every 
t > by 



E{e i{u ' Y t h) } = exp 



t( [ e^-l-i{u,y)ir(dy)) 
\J{\y\<h} J 



'{\y\<h} 

and (Nj}) is a compound Poisson process with parameters X h = ir(\y\ > h) 
and /J. h (dy) = l{| y |>h} 7r (rfy)/ 7r (|y| > h). By this remark, we want to empha- 
size that the formulation (1) is conventional and that the coefficient b in 
(1) is dependent on this choice. We will come back to this remark when 
we introduce the assumptions of the main results where we want, on the 
contrary, that they be intrinsic (see Remark 4). 

Let us recall a result about existence, uniqueness and Markovian structure 
of the solutions of (1) (see [22]). 

Theorem 1. Assume that b, a and k are locally Lipschitz functions 
with sublinear growth. Let (0,jr, (jT t ),P) be a filtered probability space sat- 
isfying the usual conditions and let Xq be a random variable on (Cl,J-,¥) 
with values in M. d . Then, for any {Tt)-Brownian motion (Wt)t>o> f or an V 
(Tt) -measurable (Zt)t>o as previously defined, the SDE (1) admits a unique 
cadlag solution (Xt)t>o with initial condition Xq- Moreover, (Xt)t>o is a 
Feller and Markov process. 

Remark 2. Levy driven SDEs are the largest subclass of SDEs driven by 
semimartingales such that the solutions have a Markovian structure. Indeed, 
a result due to Jacod and Protter (see [13]) shows that, under appropriate 
conditions on the coefficients, a stochastic process solution to a homogeneous 
SDE driven by a semimartingale is a strong Markov process if, and only if, 
the driving process is a Levy process. 

1.3. Background on approximation of invariant measures for Brownian 
diffusions. This problem has already been studied by several authors when 
{Xt) is a Brownian diffusion, that is, when k = 0. In [28] Talay approximates 
u{f) by P%{f) = l/ n ELi/(^-i)i where {X%) n denotes the Euler scheme 
with constant step 7. Denoting by the invariant distribution of the homo- 
geneous Markov chain {X%) n , he shows that 9% n ^£° v i an d that v, 
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under some uniform ellipticity and Lyapunov-type stability assumptions. (A 
Markov process (Xt) with infinitesimal generator A satisfies a Lyapunov as- 
sumption if there exists a positive function V such that V(x) — > +00 and 
limsup^4V(x) = —00 when |x| — > +00. Then, V is called a Lyapunov func- 
tion for (Xt). Under this assumption, (Xt) admits a stationary, often ergodic 
when unique, distribution. The existence of such a Lyapunov function de- 
pends on the mean-reversion of the drift and on the intensity of the diffusions 
term (see, e.g., [5, 7, 12] and [19] for literature on Lyapunov stability).) In 
this procedure, 7 and n correspond to the two types of errors that the dis- 
cretization of this long time problem generates. Practically, one cannot effi- 
ciently manage them together. Indeed, when one implements this algorithm, 
one sets a positive real 7 and then, one approximates the biased target v" 1 . 
In order to get rid of this problem, Lamberton and Pages (see [14, 15]) re- 
place the standard Euler scheme with constant step 7 with an Euler scheme 
with decreasing step 7 n . Denoting by (X n ) n >\ this Euler scheme and by 

(Vk)k>i a sequence of weights such that H n = J2k=i Vk n ^i?° +00, they de- 
fine a sequence of weighted empirical measures (v n ) and show under some 
Lyapunov assumptions (but without ellipticity assumptions) that if (r] n /^ n ) 
is nonincreasing, 

1 n 

&n(f) = jr J2vkf( x k-i) n ^±>°° v(f) p.s., 
Un k=i 

for every continuous function / with polynomial growth (see [14, 15] for 
more details and [17] for extensions). 

Remark 3. These two approaches are significantly different. Talay's 
method strongly relies on the homogeneous Markovian structure of the con- 
stant step Euler scheme and on its classical "toolbox" (irreducibility, positive 
recurrence, . . . , see, e.g., [18]). Since the Euler scheme with decreasing step 
is no longer homogeneous, Lamberton and Pages develop another method 
based on stability of Markov chains and on martingale methods which can 
be extended to a nonhomogeneous setting (see [9]). This is why they do not 
need any ellipticity assumptions on the coefficients. 

1.4. Difficulties induced by the jumps of the Levy process. In this paper 
we adapt the Lamberton and Pages approach. In order to obtain some similar 
results in the case of Levy driven SDEs, one mainly has two kinds of obstacles 
to overcome. 

From a dynamical point of view, the main difficulty comes from the mo- 
ments of the jump component. Indeed, by contrast with the case of Brownian 
motion, the jump component can have only few moments (stable processes, 
e.g.), and it then generates some instability for the SDE. 
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The second obstacle appears in the simulation of the Euler scheme. Actu- 
ally, only in some very particular cases can the jump component of a Levy 
process be simulated (compound Poisson processes, stable processes, . . .). In 
those cases, the Euler scheme [that we call exact Euler scheme and denote 
by (A)] can be built by using the true increments of (Zt). Otherwise, one 
has to study some approximate Euler schemes where we replace the incre- 
ments of Zf with some approximations that can be simulated. The canonical 
way for approximating the jump component is to truncate its small jumps. 
Let (t£ n ) n >i be a sequence of positive numbers such that u n < 1 and (u n ) 
decreases to and (Y n ) n >i be the sequence of cadlag processes defined by 

Y t n = AY ^{AY s eD n }-t [ V<dy) Vi>0 

with D n = {y £ R l , u n < \y\ < 1} and AY S = Y 8 - Y s - . The process Y n is a 
compensated compound Poisson process with parameters X n = ir(D n ) and 
fJ-n(dy) = lD n (y)^(dy)/7r(D n ). It converges locally uniformly in L 2 to Y, 
that is, for every T > 0, 

(3) e( sup \Y t - Y t n \ 2 ) n ^°° (see [4]). 

We will denote by Z n the process defined by Z n = Y n + N and by (B) 
the Euler Scheme built with its increments. The increments of Z n can be 
simulated if A n and the coefficient of the drift term can be calculated, and 
if fi n can be simulated for all n G N. This is the case for a broad class of 
Levy processes, thanks to classical techniques (rejection method, integral 
approximation,...). If there exists (u n ) n >i such that the increments of Z n 
can be simulated for all n £ N, we say that the Levy measure can be simu- 
lated. However, the simulation time of Z^ depends on the average number of 
its jumps, that is, on 7r(|y| > u n )t. When the truncation threshold tends to 
(this is necessary to approach the true increment of the jump component), 
7r(|y| > u n ) explodes as soon as the Levy measure is not finite (i.e., as soon 
as the true jump component is not a compound Poisson process). It implies 
that the simulation time of Z^ explodes for a fixed t. However, thanks to the 
decreasing step, it is possible to adapt the time step 7 n and the truncation 
threshold u n so that the expectation of the number of jumps at each time 
step remains uniformly bounded. Following the same idea, it is also possible 
to choose some steps and some truncation thresholds so that the average 
number of jumps at each time step tends to 0. In this case, approximating 
the true component by the preceding compound Poisson process stopped at 
its first jump time (the first time when it jumps) can also be efficient [see 
Scheme (C)]. 
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1.5. Construction of the procedures. Let (7 n )n>i be a decreasing se- 
quence of positive real numbers such that Hm7 n = and T n = J2k=i Ik — > 
+00 when +00. Let (U n ) n >\ be a sequence of i.i.d. square integrable 
centered Revalued random variables with Y,jj 1 = 1^. Finally, let (Z n ) n >i, 
(Zn) n >i and (Z,^) n >i be sequences of independent -valued random vari- 
ables independent of (U n ) n >i, such that, 

Zn { = ] Z, n , Z*^Zl and 2% ® Z^ Vn>l, 

with T n = inf{s > 0, |A2£| > 0}. Let x G M d . The Euler Schemes (A), (B) 
and (C) are recursively defined by Xq = Xq = Xq = x and for every n > 1, 

(A) A n+ i = A n + 7 n+ i6(A > n ) + A /7n+iO"(X n ,)C/ n+1 + K(A„)Z n+ i, 

(B) X* +1 = Af + 7 n+i6(Af ) + y7^Ia(A ? f )C/ n+ i + k(X*)Z* +1 , 

(C) A£ +1 = A^ 7 + 7n+i6(A^) + ^^o-(X%)U n+1 + /<A£)z£ +1 . 

We set T n = a(X k , k < n), = a(X^,k < n) and T% = a{X£,k < n). Let 
(%)fceN be a sequence of positive numbers such that H n = J2k=i Vk - > +00. 
For each scheme, we define a sequence of weighted empirical measures by 



(4) 



1 n 
-FT E ^x k 
nn k=l 



n 1 n 

JtY,^xb_ i and v% = r=-Y,T tk 8xa_ i . 



For a function /:IR rf i— > M, (f n (/)) can be recursively computed (so is the 
case for the two other schemes). Indeed, we have v\{f) = f(x) and for every 
n>l, 

K + l{f) = M/) + ^"(/(An+l) - Pn(/))- 

Some comments about the approximate Euler schemes. In Scheme (B), 
since (u n ) decreases to 0, we discard fewer and fewer jumps of the true 
component when n grows. We will see in Theorem 2 below that this is the 
only condition on (u n ) for the convergence of (P^)- This means that we only 
need the law of Zl l /n to be an "asymptotically good approximation" of the law 
of Zry n . Yet, as previously mentioned, there is a hidden constraint induced 
by the simulation time which is proportional to the average number vr(|y| > 
u n )ln of jumps of Z n on [0,7„]. In practice, we require (vr(|?/| > u n )^ n ) to 
be bounded. 

Furthermore, if vr(|y| > u n )^ n — > (i.e., the average number of jumps at 
each step tends to 0), we will see that the first jump of Z n is all that matters 
for the convergence of the empirical measure. This means that Scheme (C) 
becomes efficient. 
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1.6. Notations. Throughout this paper, every positive real constant is 
denoted by C (it may vary from line to line). We denote the usual scalar 
product on by (•, •) and the Euclidean norm by | • |. For any d x I real 
matrix M, we define ||M|| = sup/M^i |Mx|/|x|. For a symmetric d x d real 
matrix M, we set Am = max(0, Ai, . . . , A<j), where Ai, . . . , A^ denote the eigen- 
values of M. For every x G M d , 

(5) Mx® 2 = x*Mx<\ M \x\ 2 . 

We denote by Ct>(M. d ) [resp. Co(R rf )] the set of bounded continuous functions 
on M. d with values in M. (resp. continuous functions that go to at infinity) 
and Cf^(R rf ), the set of C 2 -functions on M. d with values in M and compact 
support. One says that / is a p-H61der function on E with values in F (where 
E and F are normed vector spaces) if 

m \\f(x)-f(y)\\ F ^ 

[fjp = SUP - |jp < +oo. 

x,y£E \\X D\\e 

Finally, we say that V : M. d i— ► Wj_ is an EQ-function (Essentially Quadratic 
function) if V is C 2 , limV(x) = +oo when \x\ — > +oo, |VV| < CVV and 
D 2 V is bounded. [In particular, V given by V(x) = p + Sx® 2 , where p is a 
positive number and S is a definite and positive symmetric dx d real matrix, 
is an EQ-function.] For p > 0, one checks that ||D 2 (V P )|| < CV p ~ l , that V p 
is a 2p-H61der function if p < 1/2, and that V p ~ l W is a (2p — 1)-H61der 
function if p G (1/2, 1) (see Lemma 3). Hence, A p and c p given by 



(6) 



A p :— — sup Xyi- PD 2{vp)(x) an d 



P 



2p 



ifpG (0,1/2], 



[V^Whp-i, if pe (1/2,1] 



are finite positive numbers. 

1.7. Organization of the paper. The main results (Theorems 2 and 3) are 
stated in Section 2 and are proved in Sections 3, 4 and 5. First, we focus 
on the proof of these theorems for the exact Euler Scheme (A): in Section 
3 we prove the almost sure tightness of (p n ) and in Section 4 we establish 
that every weak limiting distribution of (p n ) is invariant for the SDE (1). 
Second, in Section 5 we point out the main differences which arise in the 
proofs when considering the approximate Euler Schemes (B) and (C). In 
Section 6 we show that the almost sure central limit theorem for symmetric 
stable processes (see [3]) can be obtained as a consequence of our main 
theorems. Finally, in Section 7 we simulate the procedure on some concrete 
examples. 
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2. Main results. In Theorem 2 we obtain a result under simple condi- 
tions on the steps and on the weights. In Theorem 3 we show that, under 
more stringent conditions on the steps and on the weights, some assump- 
tions on the coefficients of the SDE can be relaxed. Let us introduce the 
joint assumptions. First, we state some assumptions on the moments of the 
Levy measure at +oo and 0: 

(H p ) : / *(dy)\v? p < +oc, (H 2 ) : / Tr(dy)\y\ 2 « < +oo, 
■%I>1 -%I<1 

where p is a positive real number and q £ [0, 1] . 

Assumption (Hp) is satisfied if, and only if, ¥,\Zt\ 2p < +oo for every t > 
(see [1], Theorem 6.1). By the compensation formula (see [4]), (H 2 ) is 
satisfied if and only if E{X)o<s<t l^^tl } < that is, if and only if (Y t ) 
has locally 2q-integrable variation. We recall that (H 2 ) is always satisfied for 
q = 1 since //u<h \y\ 2n (dy) < °o for any Levy measure it. 

Now, we introduce the Lyapunov assumption on the coefficients of the 
SDE and on -k denoted by (S a>P)q ). The parameter a specifies the intensity 
of the mean-reversion. We denote by b the function defined by 
(b, ifp<l/2<g, 

b-K yn(dy), Hp,q< 1/2, 

J{\v\<i} 

b + K yn(dy), ifp>l/2. 

J{\y\>i} 

The function b plays the role of the global drift of the dynamical system 
resulting from b and from the jump component (see Remark 4 for more 
precisions). Let a& (0, 1], p > and q £ [0, 1]. 

Assumption (S a ,p, q ). There exists an EQ-function V such that: 

1. Growth control: \b\ 2 < CV a , 

Tr(<7cr*) + ||k|| 2 (p v< ?) < CV a+p -\ if p < 1, 
Tr{aa*) + \\K\\ 2 <CV a , if p > 1. 

2. Mean-reversion: there exist /3 G 1, a > such that, (VV, b) +4> Piq (a, k, it, 
V) < (5 — aV a , where ^ Pi g is given by 

<j) M (a,K,iT,V) 

c p m 2 p t n\\K\\ 2p V 1 ~n q < p , ifp<l, 
Ai(Tr(a(j*)+m 2 ^|| K || 2 ), if p = l, 

dpX p [ Tr(aa*) + ™2,7r|M| 2 + e p m 2p ,n _ 1 j , if p > 1, 
with m r7r = / \y\ r ir(dy), d p = 2( 2 ( p_1 )~ 1 )+ , c p and X p given by (6), and e p = 
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Assumption (S ajPjq ).2 can be viewed as a discretized version of "AV P < 
(3 — aV a+p ~ 1 where A is the infinitesimal generator of (Xt) defined on a 
subset V(A) of C 2 (R d ) by 

A fix) = (Vf,b)(x) + \Ti(a*D 2 fa)(x) 

(7) 

+ / {f(x + K(x)y)-f(x)-(Vf(x),K(x)y)l { \ y \< l} )ir(dy). 



Furthermore, one can check that if Assumption (S a ,p,q) is fulfilled, then there 
exist £ R and a > such that U AV P < (3 - aP^" 1 ." This means that if 
V is the function whose existence is required in Assumption (S a ,p,q), then 
V p is a Lyapunov function for the stochastic process (Xt) and for the Euler 
scheme (X n ). 

The left-hand side of (S a!Piq ).2 is the sum of two antagonistic compo- 
nents: (W, b) produces the mean-reverting effect (see Example 1 for con- 
crete cases), whereas the positive function (j) pq expresses the noise induced 
by the Brownian and jump components. In particular, if the following tighter 
growth control condition holds, 



(8) \b\ 2 <CV a , 



Tr(<7 a* (x)) W = + °° (y«-(i-p)+ ( x )), 

l|K(x)|| 2|a:| = +00 o(yta(x)), 



with ?]a,p,g = (pVg) x (a + p — 1) if p < 1 and r/ a , p ,g = a if p > 1, then the 
term (j) p ^ q becomes negligible and the mean-reversion assumption becomes 

(VV,6) <(3-aV a . 



Remark 4. If we had chosen to compensate the jumps smaller than 
h > rather than h = 1, the corresponding assumption would have been 
(S>a,p,q)i where (S^ p q ) is obtained from (S a! p !q ) by replacing b with b h and b 

with b h defined by 



'b h , 




b h -K 


[ 




>{\y\<h} 


b h + K 


f 




'{\y\>h} 



yn(dy), 
yx(dy), 



iip<l/2<q, 
]£p,q< 1/2, 

if p> 1/2. 



One can check that, for every h > 0, (S a p q ) <^=^ (S a , p ,q)- This means that 
these assumptions do not depend on the choice of the truncation parameter 
h. Indeed, first, it is clear that (S apq ).l <^=> (S aiPiq ).l. Second, whenp > 1/2 



or p,q < 1/2, (S a p q ).2 <^=^> (S ajP)q ).2 because in these cases, b h = b for every 
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h > 0. This can be explained by the existence of a formulation of the SDE 
that does not depend on h. Actually, when p > 1/2, we can rewrite the SDE 
(2) by replacing b h with b h and, Z{> with i* = Z%-t J { \ y \ >h} yn(dy), that 

is, we can compensate the big jumps. Since (Z^) = (Z^°) for every h > 0, 
it follows that b h = b(=b°°) for every h > 0. There also exists an intrinsic 
formulation when p,q< 1/2 because in this case, we can replace b h with b 
and Zt with Z% = Zf* + t J^ <h yyir(dy) (now, we do not compensate any 

jumps). Since (Zf 1 ) = (Z®) for every h > 0, b =b{= b°). These formulations 
can be considered as the natural formulations of the dynamical system in 
these settings. 

When p < 1/2 < q, there is no intrinsic formulation of the SDE (even if ir 
is symmetrical). Since b h depends on h, it appears that the left-hand side of 
(£>a,p,q)-2 also depends on h. However, under the growth assumption on k, 
one can check that (VV, b h ) = (W,b) + o(V a ) and it follows that the same 
conclusion still holds in this case. 

We now state our first main result. 

Theorem 2. Let a e (0, 1], p > and q G [0, 1] such that (Hp), (H^) and 
(Sa,p,qj a^e satisfied. Suppose that E{|[/i| 2 ( pV1 )} < +oo and that the sequence 
(Vn/ln)n>i is nonincreasing . Then: 

(1) If p/2 + a — 1 > 0, the sequence (v n )n>i is almost surely tight. More- 
over, if k(x) ' ' =~ o(\x\) and Tr(o~a*) + \\K\\ 2q < CV p / 2 ~*~ a ~ 1 , then every 
weak limit of this sequence is an invariant probability for the SDE (1). In 
particular, if {Xt)t>o admits a unique invariant probability v, for every con- 
tinuous function f such that f = o{V p l 2+a ~ l ), limn^oo v n (f) = v(f) a.s. 

(2) The same result holds for (z^f ) n >i- 

(3) The same result holds for {v~()n>i under the additional condition 
(9) 7r(|y|>u n ) 7n n ^±?°0. 

We present below some examples which fulfill the conditions of Theorem 
2. In the first we suppose that the dynamical system has a radial drift term 
and a noise generated by a centered jump Levy process with a Levy measure 
close to that of a symmetric stable process. In the second we suppose that 
the SDE is only driven by a jump Levy process, but we suppose that it is 
not centered. This implies that even if the SDE has seemingly no drift term, 
a mean-reverting assumption can be still satisfied. 

Example 1. Let 4> and tp be positive, bounded and continuous functions 
on R d such that 

<f>(x) = <f>(-x) Vx G 
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6 = mm6(x) > and ib = inf ib(x) > 0. 

- Rrf — {|x|>l} 

Consider (Zj) denned as in the SDE (1) with Levy measure n given by 
n(dy) = 4>(y) /\y\ d+r Xd(dy) , where r G (0, 2). When <j> = C > 0, the increments 
of (Z() can be exactly simulated because (Zf) is a symmetric Instable 
process with order r. In the other cases, and can be simulated by 
the rejection method since the density of tt is dominated by the density of 
a Pareto's law. 

Let p G [0, 2) and 6 be a continuous function defined by b(x) = —ip(x)x/\x\ p 
We consider (X t ) solution to 

(10) dX t = b(X t -) dt + K (X t -)dZ t , 

where k is a continuous function such that ||k;(x)|| 2 < C(l + |^| 2 ) e with e < 1. 
A natural candidate for the function V is V(x) = 1 + |x| 2 . Indeed, since b = b 
[because <f>(y) = <j){— y)\, one checks that there exists (3 G M such that 

(Vy(s),6(a;)) = -2^(x)|x| 2 - p < /3 -^y(x) 1 -^ 2 . 

We set a := 1 — p/2 and 

A(r) :={(p,«) G (0,+oo) x [0,1], (Bj) and (H 2 ) hold}. 

We have A(r) = (0,r/2) x (r/2,1). By (8), for every (p,q) G A(r), (S a , Piq ) is 
satisfied if e < {jp + a — l)/q = (p — p/2)/q. Hence, (u n ) is tight if there exists 
(p, q) G A(r) such that p/2 + a— 1 > 0, that is, if p > p, and (2p — p)/{2q) > e. 
If, moreover, ||k|| 2q < C(l + \x\ 2 )P/ 2+a ~ 1 , that is, if (p - p)/(2g) < e, every 
weak limit v is invariant for the SDE (10). 

It follows that if the invariant distribution v is unique, v n ==^ v a.s. as soon 
as 2p < r and e < sup{(p — p)/{2q), (p,q) G A(r)} = 1/2 — p/r. Furthermore, 
v n (f) — > v(f) a.s. for every continuous function / satisfying f(x) < C(l + 
\x\) e with 6e [0,(r/2-p)/2). 

Example 2. Let tt be a Levy measure on M such that /|j,|<i \y\^(dy) < 
+oo, S\ y \>\ \u\ 2p7r (dy) < +oo with p > 2 and / yi:{dy) > 0. Let (Z 4 °) be a real 
Levy process with characteristic function given for every t > by 

E | e i(n,Z t 0>| =exp 

For instance, (Z^) can be a subordinator with no drift term. We assume 
that k(x) = —ip{x)x/\x\ p with p G [0,2) and ip defined as in the preceding 
example. We then consider the SDE: 

dX t = K (X t - ) dZ t ° = b(X t - ) dt + «(At- ) dZ t , 
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with b(x) = k(x) J {lyl < 1} yir(dy) and Z t = Zf - t J { ^ l} yn{dy). Since p > 

1/2, b(x) = b(x) + K(x) J { ^ >1} yir(dy) = k(x) J yn(dy). Setting V(x) = l + x 2 , 
one checks that there exists G K. such that 



We set a = 1 - p/2. Let p > 2 and <? < 1/2 such that (Hp) and (H 2 ) hold. 
First, checking that as soon as p > 0, ||k(:e)|| 2 = o(l + |x| 2 ) a when \x\ — ► +oo, 
we derive from (8) that (S aiP! q) is satisfied as soon as pG (0, 2). Second, for 
every p > 2, g < 1 and a G (0, 1), one can check that p/2 + a — 1 > and 
||k(»|| 29 < C(l + |j;| 2 ) p/2+a ~ 1 . Hence, Theorem 2 applies for every p G (0,2). 

The interest of Theorem 2 lies in the facility with which it can be put to 
use in concrete situations. For instance, in Scheme (A), we only have to take 
a sequence (7 n )n>l decreasing to 0, with infinite sum and r] n = j n . The next 
theorem (Theorem 3) requires tougher conditions on the sequences (^ n ) and 
(r] n ), but it can be applied to SDEs where the coefficients do not necessarily 
verify all conditions of Theorem 2. It broadens the class of SDEs for which 
we can find an efficient procedure for the approximation of the invariant 
measure. 

Theorem 3. Let a G (0, 1], p > and q G [0, 1] such that (Hp), (H 2 ) and 
rS a , P ,qj are satisfied. Suppose that E{|[/i| 2 ( pV1 )} < +oo. Then: 

(1) Let s G (1,2] satisfying the following additional conditions when p > 



Ifp/s + a — 1 > ; there exist some sequences (7 n )n>i <ind (?7 n )n>i such that 

ipn)n>i is almost surely tight. Moreover, if k(x) = o(\x\) andTr(aa*) + 
\\n\\ 2q < CV v l s+a ~ l , then every weak limit of this sequence is an invariant 
probability for the SDE (1). In particular, if (Xt)t>o admits a unique invari- 
ant probability v , for every continuous function f such that f = o(V p ^ s+a ~ 1 ) , 
lim^oo v n (f) = a.s. 

(2) The same result holds for (v^) n >\. 

(3) The same result holds for (^) n >i under the additional condition (9). 

Remark 5. The sequences (r? n )n>i and (7 n )n>i must verify an explicit 
condition given in Proposition 1 below (see Remark 6 for a version adapted 
to polynomial steps and weights). 




l-p/2 



1/2: 



(11) 
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In the following example, we consider the same class of SDEs as in Ex- 
ample 1 in the nonintegrable case (i.e., r < 1). One can observe that the 
mean-reversion condition and the growth condition on n and on the func- 
tions / whose the procedure converges can be relaxed. We also give some 
explicit polynomial weights and steps for which Theorem 3 applies in this 
case. 



Example 3. Let p G [0,2) and r G (0,1], let b and k be continuous 
functions defined as in Example 1. Consider (Xt) solution to the SDE (10) 
and assume that the invariant measure v is unique. For s G (1,2], denote by 
(ln,s) and (r] n ,s) some sequences of steps and weights satisfying 7„ iS = Cn~ ri , 
Vn,s = Cn~ T2 with r\ < r 2 and 

0<n<2^1--^ and r 2 < 1 or 0<ri<2^1--^ and r 2 = 1. 

Then, for these choices of steps and weights, u n ==> v a.s. as soon as sp <r 
and e G [0, 1/s — p/r) (this improves the condition: 2p < r and e G [0, 1/2 — 
p/r) of Example 1). Furthermore, P n (f) — ► v {f) a - s - f° r every continuous 
function / satisfying |/(x)| < C(l + \x\) e with 9 G [0, (r/s — p)/2) (this im- 
proves the condition: 9 G [0, (r/2 — p)/2) of Example 1). 



3. Almost sure tightness of (v n (w, dx)) n £fq. The main result of this 
section is Proposition 1. We need to introduce the function f a ^ p defined for 
all s G (1,2] by 

{s, if s > 2p, 

, , P + a -\ ^A S , if,<2p. 
p/« + (o-l)/(2(pAl)) ' V 

Assume that p + a — 1 > 0. Then, s •— > f a p( s ) is a nondecreasing function 
which satisfies fi )P (s) = s for all p > and / a , P (2) = 2. Note that f a ,p(s) > 1 
if and only if s satisfies assumption (11). 

Proposition 1. Let a G (0, 1], p > and q G (0, 1] such that (Hp), (H^) 

and (S a ,p,qJ are satisfied. Assume that E{|C/i| 2 (p v1 )} < +oo and (r) n /y n )n>i 
is nonincreasing. 



(1) T/ien, 



sup^ n (T/ p/2+a - 1 )<+oo a.s. 
n>l 



Consequently, if S + a — 1 > 0, iae sequence (v n )neN is a.s. tight. 



14 



F. PANLOUP 



(2) Let s£ (1,2) such that assumption (11) is satisfied. Assume that 
(Vn)n>i and (j n )n>l are such that 

) is nonincreasinq and [ — — — ) < +oo. 

(13) 
Then, 

su Pn>i Vniyvl s+a •*•) < +oo a.s. and i/ie sequence (^ n )neN ^s s.s. ii<?/ii 
as soon as p/s + a— 1>0. 

Remark 6. If j n = Cn~ ri and n, 
tion (13) reads 

r2 < 1 and < ri < f\ 

(14) 

r2 = 1 and < ri < ri . 

The proof of Proposition 1 is organized as follows: first, in Section 3.1 (see 
Proposition 2) we establish a fundamental recursive control of the sequence 
(V p (X n )): we show that (R a , p ): There exist n £ N, a' > 0, /3' > such that 
Vn > no, 

(15) £{V^(X„ +1 )|.F n } < V p (X n ) + Tn+i^ 1 ^)^' - «V a (X n )). 

For this step, we rely on Lemma 2 that provides a control of the moments 
of the increments of the jump component in terms of p and q. 

Second, in Section 3.3 we make use of martingale techniques in order to 
derive some consequences from (R a ,p). In Lemma 5 we establish a L p -control 
of the Euler scheme with arguments close to [15] . This control is fundamental 
for the proof of Corollary 1 where we show the following property: 

(C P)S ): There exist p £ (1,2] and a sequence (ir n ) of ^"n-measurable random 
variables such that 

( 16 ) E (-^-YmV^iXn) - 7T n -x\ P } < +oo. 

This step is used to obtain a //-martingale control (see proof of Lemma 1). 
We will see in the proof of Corollary 1 that the choice of the sequence {n n ) 
depends on p and q. In particular, even if q does not appear in (C PjS ), this 
assumption indirectly depends on this parameter. The same remark holds 
for (R ap ). In the following lemma, we show that these two steps are all what 
we have to show for the proof of Proposition 1 . 

Lemma 1. Let p>0, a G (0, 1] and s 6 (1,2] such that (H p ) ; (R a , P ) and 
(C ps ) are fulfilled. Assume that E{|[/i| 2( - pV1 )} < +oo and that (rj n /j n ) is 




= Cn 2 with r\ < ri-, then assump- 
:=2(l-—!—) or 
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nonincreasing. Then, 

(17) sup^„(F p / s+a - 1 )<+oo a.s. 

n>l 

and the sequence (P n ) n >i is a.s. tight as soon as p/s + a — 1 > 0. 

Proof. By a convexity argument (see Lemma 3 of [15]), one shows that 
(Ra,p) =^ (Ra,p) f° r an P £ (0,p]. Hence, for all s € (1,2], there exists no £ N, 
d > and (3 > such that V7c > no , 

(18) E{^/ s (X,0|JW<^ /s (*fc-i) + 7^ 

For i? > 0, set e(R) = sup {N>i?} V~ a {x) and M{R) = sup^^ V J, / a-1 (ac). 
We have 

(19) yP/ 8 " 1 ^) < £(i?)y p / s+a - 1 (a;) + Af(iJ). 

Since F(x) — > +oo when |x| — > +oo (resp. since V is bounded on every 
compact set), e(R) — > when R — > +oo [resp. M(R) is finite for every 
R > 0]. Hence, for every e > 0, there exists M £ > such that < 
£ y P /s+a-i + M ^ By getting £ = Si/ (2$), a = a/2 and = 0M E , we deduce 
that VP/^iP - aV a ) < (3 - 5iV p / s+a - 1 . Hence, we derive from (18) that 

v^-\x^) < -«<y-(*.)i*-.} + g vt> „ . 

It follows that (17) holds if 

SU P (tt E -(v p/ '(x fc _ 1 )-E{yf/'(x fc )|^ fc _ 1 })] 

n>n +lV^n fc =rlQ+1 7fc ^ 

(20) 

< +oo a.s. 

We then prove (20). We decompose the above sum as follows: 
1 n 

Hn k=n 0+ l ^ k 

1 n 

= ~ E -A^/ s (X fc ) 

1 n 



16 F. PANLOUP 

where AV p / s (X k ) = V p ' s {X k ) - F p / s (AVi). First, an Abel's transform yields 
-JT E -^V p/ %X k ) = ^(^V^(X no )-^V^(X n )) 

H n t _±~i-1 7fc tin V7n 7n / 



< ^"o yp/s(x \ 



where we used in the last inequality that (jln/ln) ls nonincreasing. Hence, 

H" -v 



since fl„ +00 and ^^W S (X„ ) n =^±?°0 a.s., 



(~ E ^AV^(X k ))<+c 



(21) sup -— ^ ^AV p ' s {X k ) ) < +00 a.s. 

n>n V fc=no 

Second, one denotes by (M n ) ng N the martingale defined by 

n 

(22) M n = J2 -fr^(V p/s (X k ) - E{V*/'(X fc )|.F fc _i}). 

fc =i ^fc 

Let /? G (1,2] and (7Tfc) be a sequence of ^-measurable random variables 
such that (16) holds. We derive from the elementary inequality \u + v\ p < 
2^ 1 (|n|^ + H p ) that 

n\V p/s {X k ) -E{^/ s (X fc )|^-i}| p } 

< cn\v p,s {x k ) - vr^n + cn\n^k-i - v^i\x k ))\^ 1 }\ p } 

<CTE{|V/'(X fc )-7r fc _i|"}, 

thanks to the Jensen inequality. Hence, (C P)S ) yields J2k>i E{| AMjt| p } < +00 
a.s. Since p > 1, it follows from Chow's theorem (see [11]) that M n n -liS? 
a.s. where is finite a.s. Then, Kronecker's lemma yields 

( 23 ) TT E ^^(XO-E^^I^^O a.s. 



fc=no+ 



! 7fc 



Hence, (20) follows from (21) and (23). Finally, since \im\ x \_> +OQ F p / s+a_1 (a;) = 
+00 when p/s + a — 1 > 0, we derive from a classical tightness criteria (see, 
e.g., [9], page 41) that {u n )n>i is a.s. tight as soon as p/s + a — 1 > 0. □ 

3.1. A recursive stability relation. 

Proposition 2. Letp>0, qe[0, 1] and a 6(0,1]. Assume (H p ), (Hq) 
and fS a ,p,qj- moreover, E{|[/i| 2 ( pV1 )} < +00, iaen (R a , p ) holds. 
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The idea of the proof of Proposition 2 is to obtain an inequality of the 
following type: 

K{V p (X n+1 ) - V p (X n )\F n } < y n+lP V p - 1 (X n )$(X n ) + R n , 

where <3? = (VV, 6) + (f) Ptq (a, k, it, V) [see (S a!Pjq ).2] and R n is asymptotically 
negligible in a sense made clear in the proof. To this end, we begin by three 
lemmas. In Lemma 2 we study the behavior of the moments of (Zt) near 
0. Then, in Lemma 3, we state some properties of the derivatives of V p in 
terms of p and in the last one (Lemma 4) we control the contribution of the 
jump component on the conditional expectation (conditioned by J- n ) of the 
increment V p (X n+1 ) - V p (X n ). 

Lemma 2. (i) Let p > such that (Hp) holds. Then, for every h > 0, 
there exists a locally bounded function iph such that 

(24) Vt>0 E{\N t h \ 2p }= f \y\ 2p ^(dy)(t + ^ h (t)t 2 ). 

J\y\>h 

(ii) Let q G [0, 1] such that (ffi) holds. Then, for every h > 0, 
e{ Y t h + 1 I y-K(dy) < t f \y\ 2 H(dy), if q < 1/2, 

I J\y\<h J J\y\<h 

E{\Y t h \ 2 «}<C q t ( \y\ 2 ^(dy), ^gG (1/2,1]. 

J\y\<h 

(hi) Let p G [l,+oo) such that (Hp) holds. Then, there exists rj > 1 suc/i 
i/icrf, /or every T > 0, for every e > 0, i/iere exists C £) T,p > suc/t t/iai, 

VtG [0,T] E{|Z t | 2 f}<^| |y| 2p 7r(dy)+e) +Q, r , P t r ', 

where (Zt) is the compensated jumps process defined by Zt = Zt — t J\ y \ >x yi:(dy). 
Ln particular, E|^| 2 = t f \y\ 2 ir(dy). 

Remark 7. In this lemma we obtain, in particular, a control of the ex- 
pansion of 1 1 — ^ E{|Df| r } in the neighborhood of (where D denotes one of 
the above jump components and r, a positive number). We have the fol- 
lowing type of inequality: E{|Z)^| r } < c r t + 0(t ri ), where c r is a nonnegative 
real constant and rj > 1. In the first and in the last inequality, we minimize 
this value because it has a direct impact on the coefficients of the func- 
tion Pj(? and then, on the mean-reverting assumption (see Lemma 4 for 
details). Note that we cannot have c r = in the inequalities of Lemma 2. 
Indeed, according to the Kolmogorov criterion, a Levy process D that satis- 
fies E{|Df| r } < Ct 11 in the neighborhood of is pathwise continuous [for the 
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Brownian motion, c r = as soon as r > 2 since E{|Wi| r } = o(t r ^ 2 )]. When 
p > 1, this feature generates a specific contribution of the jump component 
on the mean-reverting assumption ((S a ,p,q).2). This contribution appears in 
<j) Pt q where there is an additional term of order 2p coming only from the 
jump component. 



Proof of Lemma 2. (i) (iV^) t >o is a compound Poisson process with 
parameters A^ = ir(\y\ > h) and /J, h (dy) = lny] > t l \Tr(dy)/ir(\y\ > h). Hence, 
(iV/ 1 ) can be written as follows: Nf 1 = 2~2 n >i^"n,^T n <t, where (R n )n>i is a 
sequence of i.i.d. r.v. with law [i h and (T n ) ng N is the sequence of the jump 
times of a Poisson process with parameter A^ independent of (Rn)n>i- We 
have 



E{\N( l \ 2p } = J2 E< 



n>l 



5> 



IV. 



where F x {t) = e~ xt ^ 



n>0 



X h m{\Rt\ 2 P}F Xh (t) 
E{\Ri\ 2 p} (n + 1)!' 



By the elementary inequality (this inequality will be usually needed in the 
sequel for the control of the moments of some sums of jumps) 



(25) 



Vai,...,a n GR%Va>0 



n 



i=l 



used with a = 2p, we obtain 



E{| Er=i ^l 2p l < + l) (2p " 1)+ E& 1 l^l 2p } _ (n + l)( 2p - x )+ 



(n + l)!E{|i?!| 2 P} 



(n + l)!E{|^i| 2 P} 



It follows that F\ h is an analytic function on M such that F\ h (0) = 1. There- 
fore, 

F Xh (t) = l + ti> h (t) with \Mt)\ <C(p,h,X h ) Wte[0,T]. 

Since E{|i?i| 2p } = ^- J{|y|>/i} \y\ 2p7T (dy), the first equality is obvious. 

(ii) If j\y\<h \y\ 2qiT {dy) < +o° w hh q < 1/2, then has locally bounded 
variations and Y t h + t J^ <h yTr(dy) =Eo<s<t^(' i - Inequality (25) with a = 
2q and the compensation formula yield 

e{ Y t h + 1 [ yTT(dy) ""I <e\y, \&Y t h \ 2 A 

= t f \y\ 2q ir(dy). 

Jy<h 
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Now, let q G (1/2, 1]. As Y h is a martingale, we derive from the Burkholder- 
Davis-Gundy (BDG) inequality (see [4]) that 

E{|y/fn<cye{( £ \ay^ j. 

The second inequality follows from inequality (25) with a = q and from the 
compensation formula. 

(iii) One first considers case p = 1. The process (Mf ) defined by Mt = 
\Zt\ 2 — t J \y\ 2 n(dy) is a martingale. Then, in particular, E{|Z^| 2 } = tj \y\ 2 ir(dy). 
Suppose now that p > 1. In order to simplify the notation, we assume that 
T < 1. The BDG inequality yields 

(26) E{|y/f" } < c p e{ ( 2 l Ai ?l 2 ) ? }- 

For every integer k > 1, M tjfe := Eo<s<i l A ^)| 2fc ~ */{| v |<fc} If l 2 ** - ^!/) is 
a martingale. By inequality (25) and the BDG inequality applied to (M^), 
we obtain 



\0<s<t 



< cU{\M t ^ 2k - 1 } + (t [ \yfir(dy)) P/2 

\ \ J{\y\<h\ J 



{\y\<h} 



<C\ 



ElA^r 1 ) M 2fc ^)) P/2 
V I Vo^<* / J V h\y\<h} J 

Set A:o = inf{/c > 1, 2 k > p}. Iterating the preceding relation yields 

e{ ( E l Ay if V } ^ CE { ( E |A^i 2fc0+1 ) P/ ) 

I \0<s<t J J I \0<s<t / J 



fc / /. \ p/2 



)t-i 



{|y|<M 



~'ir(dy)j 



By construction, p/2 k ° < 1. We then derive from inequality (25) with a 
p/2 k °, from the compensation formula and from (26) that 

C / \ p> 

E{|y/f p } < 



(27) 

< C p t / |y| 2p 7r(dy) + C^" 1 

J\\v\<h\ 
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with 771 = p/2 k °~ 1 > 1. We now consider (Zt). For every h G (0, +00), we 
have Z t = Y t h + iV/ 1 , where Nj^ = Nf — t f< y < >h yir(dy). Using the elementary 
inequality, 

(28) Vu,?;G]R + ,Va> 1 (u + v) a < u a + a2 a - l (u a ~ l v + v a ), 
we derive from (24) that 

(29) Va > 1 E{\N t h \ a } < t f \y\ a 7r(dy) + C a>h t aA2 . 

J\y\>h 

Using (28) and the independence between (Nf 1 ) and (Y t h ) also yields 
E{|Z t | 2p }<E{|iV/ l | 2p } + C'(E{|iV/ t |} 2 ^ 1 E{|y^i}+E{|y t /l | 2p }). 
Since EjlY/ 1 ! 2 } = t J{\ y \<h} \y\ 27r (dy), we derive from the Jensen inequality 

that E{|y^|} < C h y/i. Hence, by (27) and (29), it follows that, for every 
h > and t < T, 

E{\Z t \ 2p }<t I \y\ 2 V7T(dy) + Cl h t 3 ^+C 2 p t [ \y\ 2p 7r(dy) 
J{\y\>h} " P h\y\<h} 

with 771 > 1, Cp h > and C 2 > 0. Let e be a positive number. As C 2 does not 
depend on h, and Jiyi<^ \y\ 2p ir(dy) — ► when h — > 0, we can choose h e > 



such that C 2 L y \ <h£ \y\ 2p ir{dy) < e. That yields the announced inequality 



□ 

Lemma 3. Let V be an EQ-function defined on R d . Then: 

(a) If p £ [0, 1/2], U p is a-Holder for any a £ [2p, 1] and if p £ (0, 1], 
V(U p ) is a-Holder for any a € [2p - 1, 1] H (0, 1] . 

(b) Let x, y £ R d and £ G + y] and sei u = min{U(x), 3; G // 
P<1, 

(30) ^ 2 (n(£)y® 2 <P^%M 2 - 
//, moreover, \y\ < (1 — e) wzi/i e G (0, 1], i/ien, 

(31) lD 2 (VV)(0y® 2 <p\e 2 ^VP-\x)\y\ 2 . 
Ifp>l, 

(32) ii9 2 (n(Oy 02 <pA p 2( 2 ^ 1 )- 1 )+(U^ 1 (x) + [v^]Jy| 2 ( p - 1 ))|n| 2 . 

Proof. Consider a continuous function / :M. d 1— ► R. Let a G (0,1] such 
that I /I 1 /" is Lipschitz. Then, / is an a-Holder function. This argument 
yields (a) (see [21] for details). Now, let us pass to (b). We have 

(33) D 2 (V P ) = P V^ 1 [D 2 V + (p - 1)^|- W ' 
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where (W ® V7)ij = (W^W).,-. Since V^" 1 < v^ 1 if p < 1, we de- 
rive (30) from relations (5) and (6). For (31), we consider £ = a; + 9y with 

G [0, 1] and \y\ < (1 — e) ^ ■ As is a Lipschitz function, 

W(0 > \/V(x) - [W^lvl > eW(x) y p_1 (0 < e^" 1 ^" 1 ^)- 

Hence, inequality (31) follows from (6). If p> 1, 

V^(0<^(a:) + [V^1ill/| =>> ^ p " 1 (e)<(v / F(x) + [v / T 7 ] 1 |y|) 2(p " 1) . 

We then derive (32) from (25) [with a = 2{p — 1) and n = 2] and from (6). 

□ 

Lemma 4. Let p G (0,1), gG [0,1] and aG (0,1]. Assume (Hp), (H 2 ) 
and ^Sa iPj qj.i. Then, /or every e > 0, there exists h £ G [0,+oo], T £ > and 
C e > suc/t t/iat /or every x, z G /or every t <T £ , 

¥.{V p (z + k(x)Z^)-VP(z)} 

(34) 

<i(pc p | |y|*7r(dy)l { ,< p} ||K(x)||* + eV^ a - 1 (x)+C e ), 
mi/j c p given oy (6), h £ G (0, 1] if p < 1/2 < q, h £ = if p,q < 1/2 and 

h £ = +oo if pe (1/2,1). 

PROOF. Set A(z,x, C7) = yp(z + k(x)C/) - ^(z). We first consider the 
case p < 1/2 and g > 1/2. Let /i G (0,oo). Since Zj/ = Y t h + Af, we can 
decompose A(z, follows: 

A(z, x, Z t h ) = A(z + n(x)N t h , x, Y t h ) + A(z, x, N t h ). 

One controls each term of the right-hand side. On the one hand, as V p is 
2p-H61der with constant [y p ]2p =pc p [see (6)], we deduce from Lemma 2(i) 
that 

E{A(z,x,N t h )}<pc p \\K(x)\\ 2 PE{\N t h \ 2 v} 

(35) 

<pc p \y\ 2 PTr(dy)Mx)\\ 2 v(t + Mt)t 2 ), 

J\y\>h 

where iph is a locally bounded function. On the other hand, we set z, = 
z + k(x)NJ/. By the Taylor formula, 

A(z,x,Y t h ) = (V(VV)(z), K (x)Y t h ) + (V(V)(0 - V(V)(5),«(x)y t h > 

with £ G [5, 5 + K(x)y t h ]. As (iV/ 1 ) and (y/ 1 ) are independent and Y t h is cen- 
tered, E{{V(V p )(z),K(x)Y t h )} = 0. By Lemma 3, V^W = V(y p )/p is 
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(2q - 1)-H61der (because 2q — 1 £ [2p — 1,1] D (0, 1] in this case). Then, it 
follows from Lemma 2(h). 2 that 

E{A(z + K (x)7v/ i ,x,y^)}<p[^ 1 vy] 2g _i|| K (x)|| 2( ?E{|y^| 2 n 

(36) 

<C\\K{x)fn / |y| 2 %(dy). 

J\y\<h 

Let e > 0. First, by (S a ,p, q ).l, || 2 ^ < CV p+a ~ 1 . Then, using that 

f\v\<h \y\ 2q7r (dy) ~^ when h — ► 0, we can hx /i e G (0, 1] such that 

(37) E{A(z + K(x)Nt,x,Y t h *)} < ^W^-\x). 

Second, since tph e is locally bounded, it follows from (35) that there exists 
C\ such that, for every t < 1, E{A(z, x, E )} < C^t\\K(x)\\ 2p . Now, as p < q, 
for every 5 > 0, there exists Cf > such that ||k(x)|| 2p < ( 5V p+a ~ 1 + Cf [see 
(19) for similar arguments]. Hence, setting 6 e = ej{2C\) yields 

(38) E{A(z,x,Nh} <t(^VP+ a -\x) + C e ) 

with C e = CgC| . Then, adding up (37) and (38) yields the result when p < 
1/2 <q. 

When p, g < 1/2, we deal with (Z t ) = (Z°). For every h>0,Z t = Y t h +N£>, 
where Y t h = Y t h + t fn y \<h} V^idy). Hence, for every h > 0, 

A(z, x, Z t ) = A(z + K(x)N t h ,x, Y t h ) + A(z, x, iV t ft ). 

If g < p, tt satisfies (H 2 ). Since p< 1/2, V p is 2p-H61der. Therefore, by 
Lemma 2(h). 1, 

E{A(z + K(x)N t \x,Y t h )} <pc p tMx)\\ 2p [ \y\ 2p 7r(dy). 

J\y\<h 

By summing up this inequality and (35), we deduce (34). When p <q < 1/2, 
we use that V p is 2q- Holder (see Lemma 3) and a proof analogous to the 
case p < 1/2 < q yields the result. 

Finally, we consider the case p > 1/2 where we deal with Zt = Z^ . For 
every h > 0, we have Z t = Y t h + Nj 1 , where Nj* = Nj 1 - fn y>h x y^{dy). For 

every h> 0, A(z,x,Zt) can be written as follows: 

(39) A(z, x, Z t ) = A{z + k(x)N[\x, Y t h ) + A{z, x, N?). 

One the one hand, by the same process as that used for (36) and by inequality 
(29), we have 

E{A(z,x,iV t h )}<p[V^ 1 VV]^i||«(x)||*E{|^|*} 

<t\\K(x)\\ 2p ( f \y\ 2p K{dy) + C h t 2p - 1 ). 

\J\y\>h J 
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On the other hand, by Lemma 3, V P ~ 1 X7V is (2(p V q) — 1)-H61der. Hence, 
(36) is still valid in this case if we replace q with pV q. By using this control 
for the first term of the right-hand side of (39), we obtain if q < p 

K{A(z,x,Z t )} 

<t\\ K (x)\\ 2p ( f \y\ 2p ir(dy)+ [ \y\ 2 ^{dy) + C h t 2p - 1 ). 

\J{\y\>h} J\y\<h J 

The result follows in this case. When q> p, the sequel of the proof is similar 
to the case p < 1/2 < q. □ 



3.2. Proof of Proposition 2. For this proof, one needs to study separately 
the p < 1 and p > 1 cases. We detail the first case. When p > 1, we briefly 
indicate the process of the proof which is close to that of Lemma 3 of [15]. 

Case p < 1. For h > 1, we set Z% = Z n - r YnJ{i < \ y \<h}y'K{dy) and for 

he (0,1), 2^ = Z n + -i n f {h< \ y \< l} yir(dy). liq< 1/2 (resp. p> 1/2), we can 
take h = (resp. h = +oo). Thus, we can write 

3 

AX n+1 := X n+1 -X n = E A ^n+i,fc with AX£u,i = ln+ib h (X n ), 

(40) 

AX^ 2 = ^^a(X n )U n+1 and AX* +lj3 = «(X n )Z* +1 . 

The idea is to study the difference V p (X n+ i) — V p (X n ) as the sum of 
three terms that correspond to the above decomposition. For k = 1,2,3, 
set X^ +l k = X% + J2i=i ^ X n+i,i- 

(i) First term: There exists n\ G N such that, for every n > n± 

nv p {x h n+ i,i)-w(x n )\T n } 

(41) h 
Indeed, from Taylor's formula, 

(X7V h h \ 1 

V p {X h n+hl ) - V p (X n ) =p ln+1 ^2_±( Xn ) + -D 2 (V p )(e n+1 )(AX^ 1A r 2 , 

where G [X n , X n + j n+1 b h (X n )]. Set x = X n and y = j n+ ib h (X n ). Since 
7n and \b h \ < CW by (S a , P ,q).l, there exists n± G N such that, 

for n>ni, \y\ < ^f^ffi a.s. Thus, we can apply the second inequality of 

2[vv]i 

Lemma 3(b) with e = 1/2 and deduce (41) from (S aiPjq ).l. 

(ii) Second term: For every e > 0, there exists ii2, £ G N such that, for 
every n > n2, £ , 

(42) E{V p (X* +h2 ) - V*>{X*)\T n } < e ln ^V a ^-\X n ) + C l eln+l . 
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Let us prove this inequality. Since E{£/ n +i|.F n } = 0, we deduce from Taylor's 
formula that 

nv p (xn +1>2 ) - n^+i,i)i^4 = \nD 2 {v p ){e n+ i){^x n+l ^\T n } 

with £ 2 +1 G [X% +ltl ;X% +lj2 \. Set x = X% +1>1 and y = s /j^a{x)U n j rl . By 
(S a ,p,q)-1 5 11^(^)11 < C a y/V{x) because p + a — 1 < 1. Then, the conditions 
of (31) are satisfied with e = 1/2 if |f/ n+ i| < p n +i = 1 / (2C<j [vr] i y^n+i ) • 
Therefore, 

E{ J D 2 (^)(^ +1 )(Ax n+1 , 2 )^i {|f/ri+ll < Pn+l} |j:4 
<c 7n+1 yf- 1 (X)^(^*)(^n) 

<C ln+1 V a+2 ^' l \X n ) 
since Tr(cr<7*) < Cy p+a " 1 when p < 1. By (30) and (S a , p , q ).l, we also have 

E{D 2 (F)(^ +1 )(AX n+ i,2) 02 i { |(/„ +1 |> Pn+l} |^} < a+nn+i^- 1 ^), 

where S n = ~E{\U n \ 2 l{\u n \ >Pn }}. Now, let e > 0. First, since a + 2(p— 1) < a + 
p- 1 when p < 1, there exists C £ > such that V^' 1 Tr(cro-*) < e^+P" 1 + C e 
[see (19) for similar arguments]. Second, since p n — > +oo, 5 n — > 0. Thus, there 
exists n 2 , e G N such that, for every n > n 2 , e , C5 n V p+a - 1 < eV p+a ~ 1 . The 
combination of these two arguments yields (42). 

(iii) Third term: For every e > 0, there exists h £ G [0,oo], Cf > and n%^ £ 
such that, for all n > n^^, 

, , E{^(^ lj3 )-^(X^ 1>2 )|^} 
(43) 

< ln+1 (pc p J |y| 2 ^(dy)l {g < p} || K (X)|| 2p + ^ +a - 1 (X n ) + C e 2 j 

with /i e G (0, 1] if p < 1/2 < q, h e = if p,g < 1/2 and /i £ = +oo if p G 
(1/2,1). This step is a consequence of Lemma 4: since U n+ \ and are 
independent, we have 

E{VV(Xt+l, 3 ) ~ V P (XHU,2)\Fn} =M{Gh e (Xt+l,2,Xn)\Fnh 

where G h (z, x) = E{V p (z + k(x)Z?) - V p (z)}. Then, Lemma 4 yields (43). 

We can now prove the proposition. Let e > 0. By adding (41), (42) and 
(43) and using that 7 2 < e^ n for sufficiently large n (since 7 n — > 0), we obtain 
that there exists n e G N, h e > and C e > such that, for every n>n e , 

E{V p (X n+1 )\F n } 

< V p {X n ) + ln+l {eV p+a -\X n ) + C £ ) 

(44) 

+ ln+l pV p -\X n ) 



[{VV,b hs ) + l q < P c p J lyl^Tr^llKll 2 ^ 1 "^^). 
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When p,q < 1/2 (resp. p > 1/2), h £ = (resp. h £ = oo). We deduce that 
(VV,b h -) = (W,b) because 6 = 6° (resp. b = b°°) when p,g < 1/2 
(resp. p > 1/2). We then recognize the left-hand side of (S aiPjq ).2 in (44). 
When p < 1/2 < q, h £ G (0, oo) and (VV, b h ^) = (VV, 6) + $ft £ , where = 
(W(x),k(x) f {hE<M < l} yir(dy)). Therefore, by (S a) p,q)-2, we obtain 

E{V p (X n+1 )\F n } 

< V p (X n )+ ln+lP V p ~ l (X n )(p - aV a (X n )) 

+ 7n+i (eV p+a ~ 1 (X n ) + C £ + l {p < 1/2<q}P V p -\X n )<$> hE (X n )). 

When p,q< 1/2 or p > 1/2, we set e = pa/2 and obtain (R a , P ) with /?' = 
p(5 + C e /v p ~ 1 and a! = pa/2. When p < 1/2 < q, by (S a , p , q ).l, one checks 
that, for every e > 0, there exists C £ > such that V p ~ l \$ he \ < e y p+a_1 + C £ 
and the result follows. 

Case p > 1. Thanks to Taylor's formula, 

y p (x n+1 ) = v p (x n ) + ln+1 (v(v p )(x n ),Ax n+1 ) 

+ ±D 2 (V p )(Z n+l )(AX n+1 )® 2 , 

where £ n+ i G [X n ,X n +i] and 

AX n+ i = A" n+ i - X n = 7 n+ i6°°(X n ) + y/j n+ icr(X n )U n+ i + ^(X^Z^ 

with = Z n — 7n//|j,i>n yn(dy)- Using that b = b°° in this case and that 
E{C/ n+ i|JP n } = E{Z~ +1 |^ n } = yields 

E{y p (x n+ i)|^n} = +P7n+i^ p - 1 (x n ,)(vy,6)(x ri ) 

+ iE{D 2 (^)(e„ + i)(Ax n+1 )^ 2 i^„}. 

The sequel of the proof consists in studying the last term of this equality. 
The main tools for this are the last inequality of Lemma 3 which provides a 
control of D 2 {V p )(£ n+ i){AX n+ i)® 2 and Lemma 2(iii), which gives a control 
of the moments of the jump component (see [21] for details or [15] for a 
similar proof). 

3.3. Consequences of Proposition 2. In Proposition 2, we established 
(R a , p ). According to Lemma 1, it suffices now to prove (C PjS ). This property 
is established in Corollary 1 and is a consequence of Proposition 2 [under 
additional assumptions on (p/ n ) and (%) when s < 2]. More precisely, we 
first show in Lemma 5 that a supermartingale property can be derived from 
(R a , p ) and that this property provides an L p+a ~ 1 -control of the sequence 
(V(X n )) [see (45)]. Second, we show in Corollary 1 that we can derive (C PjS ) 
from this lemma. 
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Lemma 5. Let a G (0, 1] andp > 0. Assume (Hp) and (R a)P ). Let {d n ) n g® 
be a nonincreasing sequence of nonnegative numbers such that Xm>l ®nln < 
Co. Then, there exists no > 0, a > and /3 > snc/i that (S n ) n > no defined 
by 

n 

s n = e n v p (x n ) + e klk v p+a -\x k ^) + $ ]T o klk 

k=l k>n 

is a nonnegative L 1 -supermartingale. In particular, 



^0 n7 „E{y^ 1 (^n-i)}<+oo and K{V P (X n )} n ^°° O 




(45) 

Proof. Since b, a and k have sublinear growth and Z n G L 2p for every 
n > 1, we can check by induction that, for every n > 0, V p (X n ) is integrable. 
Denote by (A n ) n >i the sequence of martingale increments defined by A n — 
V p (X n ) - E{yP(X n )|^" n _i}. By (R a , P ), there exists n G N such that, for 
every n > no, 

6n+lV P (X n+ i) 

< 9 n+1 A n+1 + 6 n+1 E{V p (X n+1 )\F n } 

< 9 n+1 A n+1 + 6 n+l (V p (X n ) + ln+l V p - 1 (X n ){[3' - a'V a (X n ))). 

By the same argument as in (19), one can find a > and (3 > such that 
V p ~ l (f3' — a'V a ) <$ — aV p+a ~ 1 . Since (6 n ) is nonincreasing, we deduce that 

9 n+1 (V p (X n+1 ) + a ln+1 V p+a -\X n )) 

< e n v p {x n ) + e n+1 A n+1 + o n+lln+1 p. 

Adding Ci &Ek=i^kV p+a ' l (X k ^) + /3£ fe>n+1 hlk to both sides of the 
inequality yields 

S'n+l < S n + n+ i/\ n+ i =>• E{S n+ i | Fn} <S n Vn > n . 

Since S no G L 1 , it follows that (S n ) n > no is a nonnegative supermartingale 
and then, that supE{5 n } < +oo. The result is obvious. □ 

Corollary 1. Let aG (0,1], p>0 and gG [0,1]. Assume (Hp), (H 2 ) 
and (S a .p,q,). If lE{|f/i| 2( - pV1 - > } < +co and (ry n /7 n ,) n6 N is nonincreasing, 

( 46 ) E f #^) 2 E{|^/ 2 (X„) - ^/ 2 (X n _! + 7n 6(A > n _ 1 ))| 2 } < +oo. 
In particular, (C Pi 2) ZioWs wii/i p = 2 and ix n = V p l 2 {X n + ^ n b(X n )) . 
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Furthermore, if conditions (11) and (13) are satisfied for s G (1,2), 

E (-^y^niV^iXn) - VPl S {X n _ X +7nB(^n-l))| /aJ,W } 

(47) " 

< +oo. 

In particular, (C p>s ) ZioWs i/ftfft /) = f a ,p(s) and ir n = V p f s (X n + 7 n 6(X n )). 

Proof. Let us begin the proof by two useful remarks. First, (46) is a 
particular case of (47) since / a , P (2) = 2 and (13) is always satisfied in this 
case. Indeed, as (rin/lnlnG^ is nonincreasing, so is (■£-( „ Vn ~— ) 2 ), and 

In tl n y /n 



with Aff n = i/ n — H n —\. Then, it suffices to prove (47). Second, by Lemma 
5 applied with 9 n = ^( H ^" )f a ^ s \ we have 



(49) (jT ! ^) faAs) ny p+a ^(x n - 1 )} < +oo. 

„>1 V^nV7n/ 

Hence, one checks that (47) holds as soon as 

E{\V*'"(X n ) - ^(AVx + 7n &(A > n _ 1 ))| / <-M} 

(50) 

< c 7i { Q ^ (s)/2 E{y p+a - 1 (A > ri _ 1 )}. 

Thus, we only need to prove (50). We inspect the p/s < 1/2 and p/s > 1/2 
cases successively. 

Case p/s < 1/2. In this case, f a ,p(s) = s. We keep the notation intro- 
duced in (40), with h = 1 if p < 1/2 < q, h = if p, q < 1/2 and /i = +oo if 
p > 1/2, and derive from (25) that 

\V p/s (X n ) - V^ s (X n ^ + 7n 6(X ri _ 1 ))| s 

(51) 

< c\v p /%x^ 2 ) - vp/'ix^i* + c|y p / s (x n ) - v p ' s {x^ 2 )\ s . 

We study successively the two right-hand side members. First, by the Taylor 
formula, 

\v p/s (x^ 2 ) - v p i s {x^)\ s < cri /2 \(v p ' s - 1 vv(£)MXn-i)u n )\' 

with £ [X^i;X^ 2 ]- The function V p / s ~ l W is bounded. Hence, since 
||o"|| s < CTr(o"(7*) (because s < 2), we derive from (S a . p>q ).l that 

(52) E{\V p / s (X*) - V^iX*)]^!} < C^V^-'iX^). 
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Second, since U n and Z% are independent, 

n\v p/s {x n ) - V^ s (X^ 2 )\ s } = E{T,(A^ 2 ,X n -i,7n)|^-i} 

where T h (z,x,j) = E{\V p / s (z + k(x)Z!}) - V p / S {z)\ s }. 

By (51) and (52), one checks that (50) holds if there exists C > such that, 
for every z, x E M. d and 7 < 71, 

(53) T^x^KC^V^-^x), 

where h = 1 (resp. h = 0, resp. h = +00) if p < 1/2 < q (resp. if p, q < 1/2, 
resp. if p > 1/2). Then, it suffices to prove (53). First, when p < 1/2 < q, 
we have = Z 7 = + iV 7 . On the one hand, since V p l s is a 2p/s-Hdlder 
function [see Lemma 3(a)], it follows from (S aiPiq ).l and Lemma 2(i) that 

^{\V p /%z + k{x)N 1 )-V p /%z)\ s }<C\\k{x)\\ 2p ¥,{\N 1 \ 2p } 

(54) 

<C-fV p+a ~ l (x). 

On the other hand, V p l s is a 2g/s-H61der function when q/s < 1/2 (be- 
cause 2p/s < 2q/s < 1 in this case). Hence, using this property if q/s < 1/2 
and the Taylor formula if q/s > 1/2 yields 



E{|y p / S (z + k(x)(JV~ + Yy)) - V p/s (z + K(x)iV 7 )| s } 

(55) 

"|/c(x)|| 2 9E{|y 7 | 2 9}, ifg/s<l/2, 
EIK^-Wfe),^)^)! 5 }, if g/ S > 1/2, 

with £ 2 £ [z+/e(x)JV 7 ,s+«(x)(JV 7 + l^)]. By Lemma 2(ii).2, E{|Y 7 | 2( >} < C7. 
It follows from Jensen's inequality that 

(56) E{\V p / s (z + k(x)Z 7 

) _ V p/s (z + k(x)N 7 )\ s } < C , 7 s/(29)a1 ||k(x)|| sA2( '. 

One checks that ||/c|| sA2 « < CV p+a ~ l under (S a , p , q ).l and that 7 + 7 s /( 2 ^ 1 < 
C7 S//2 for every 7 < 71. Hence, summing up (54) and (56) and using (25) 
yields (53) (with h = l) when p < 1/2 < q. 

When p,q < 1/2 (resp. p > 1/2), we have to check that (53) holds with 
h = (resp. with h = +00). Then, we need to use a decomposition of the 
jump component adapted to the value of h. We split up Z® (resp. Z°°) as 

follows: Z° = Y 7 + 7V 7 with Y 7 = Y 7 + 7 J {lyl < 1} yTr(dy) [resp. Z™ = Y 7 + 7V 7 
with 7V 7 = iV 7 — 7 J/u>i} 2/ 7r (^2/)]- Then, when p, q < 1/2 (resp. p > 1/2), the 
idea is to replace Y 7 with Yy (resp. 7V 7 with 7V 7 ) in the left-hand sides of (54) 
and (55) and to derive some adapted controls from Lemma 2 and inequality 
(29). Since the proof is close to that of the p < 1/2 < q case, we leave it to 
the reader. 
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Case p/s > 1/2. Since p > 1/2, we use the notation introduced in (40) 
with h = +oo. We recall that X n _i + 7 n 6(A > n _i) = X^. Then, applying the 
following inequality, 

(57) Vu,w>0,Vq>1 \u a -v a \ <C a {\u-v\u a ~ l + \u-v\ a ) 



with u = y/V(X n ), v = y/V(X^ 1 ) and a = (2p)/s, we obtain 

\v p ' s {x n ) - v p '\x™ x )\ < c\W(x n ) - Wix^v^-^ix^x) 

+ c\W(x n )-Vv(x™ 1 )\ 2p/s . 

We deduce from (S ajPiq ).l that 

IX _xooi < f CVia+P ~ m2P KXr^l)(V^\Un\ + \Z^\), ifp<l, 

i „ ^-\cv°/ 2 (x n -i)(y/%\u n \ + \zz>\), if P >l. 

Since y 7 ^ is Lipschitz, one then checks that 
\V*' 8 (X n )-V*' 8 (X% 1 )\ 

< CV r (X n ^)(^\U n \ + \Z™\ + Yn /s \U n \ 2p/s + |^| 2p/s ), 

'a + p — 1 



where r ■ 




if p < 1, 
if p>\. 



One derives from Lemma 2 (hi) and from the Jensen inequality that, for 
a > 0, E{|Z£°| a } = 0(7i a/2)A1 ). Therefore, since 2p/s > 1/2 and / 0)P (s) < 2, 
we have 

^{(^1^1 + +^ /S \U n \ 2 P/ S + |Z~|2p/')/».p«} = 0( 7 £*W/ 2 ). 

Second, one deduces from the definition of f a ,p that rf a , P (s) < a + p— 1. 
Therefore, inequality (50) follows. □ 

By Lemma 1, Corollary 1 concludes the proof of Proposition 1 and then, 
the part which is concerned with the tightness of {9 n {uj,dx)) n >i. The only 
thing left to prove the theorem for Scheme (A) is thus to identify the limit. 
This is the aim of the next section. 



4. Identification of the weak limits of (u n (u), dx)) n > x . In this section 
we show that every weak limiting distribution of (i> n {w, dx)) n >i is invariant 
for (Xt)t>o- For this purpose, we will rely on the Echeverria- Weiss theorem 
(see [7], page 238, [14] and [16]). This is a criterion for invariance based 
on the infinitesimal generator A of (X t ) defined by (7). By the Echeverria- 
Weiss theorem, we know that if j4(C^(IR rf )) C Co(M d ), a probability v is 
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invariant for the SDE if for every / G C]^(M. d ), v{Af) = 0. One can check that 
A(C^(R d )) cC (M d ) if ||«(a;)|| =o(|x|) when \x\ ->■ +00 (and this condition 
cannot be improved in general). Hence, under this condition on k, it follows 
that every weak limiting distribution of {v n ) is invariant if for every / G 
C^(K d ), v n (Af) — ► 0. The main result of this section is then the following 
proposition. 



Proposition 3. Let a G (0, 1], p > and q G [0, 1]. Assume (Hp), (H 2 ), 

(S a ,p,q)-1- Assume that \\k(x)\\ = o(\x\) and that (rj n /j n ) n >i is nonin- 
creasing. If, moreover, 

supP n (||K|| 2lJ + Tr(c7(J*)) < 00 and 

n>l 

(58) 

^i E{ ya+p-l ( X jfe _ 1 )}<+00 ! 
fc >l n k~ik 

then, 

(59) V/GC|(M d ), a.s., J Af dv n n ^? 0. 

Consequently, a.s., every weak limiting distribution of [v n (uj, dx)) n >i is in- 
variant for the SDE (1). 



Remark 8. This proposition is sufficient to conclude the proof because 
the two assumptions in (58) hold under the assumptions of Theorem 2 (resp. 
Theorem 3). Indeed, since ||k|| 2<? + Tr(aa*) < CV v l sJra - x with s = 2 in The- 
orem 2 [resp. with s satisfying (11) in Theorem 3], the first is a consequence 
of Proposition 1. Likewise, the second is a consequence of Lemma 5 applied 
with 9 n = (r] n /(H nln )) 2 [see (48)]. 



4.1. Proof of Proposition 3. The proof of Proposition 3 is built in two 
successive steps that are represented by Propositions 4 and 5. In Proposi- 
tion 4 we claim that showing that v n (Af) — > a.s. is equivalent to showing 
that l/ff ft £JU(W7*)E{/(X fc ) - /(X fc _i)|^ife_i} a.s. Then, in Propo- 
sition 5 we show that this last term does tend to 0. 



Proposition 4. Assume that the assumptions of Proposition 3 are ful- 
filled. Then, for every f G C 2 K (R d ), 



(60) 



a.s. 
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We begin the proof by a technical lemma. 

Lemma 6. Let <J> :M. d i— > M! be a continuous function with compact sup- 
port, ^ :M. d ^R +) a locally bounded function, (/i]%e[o,i] an d (^2)fe[o,i] two 
families of Borel functions defined on M. d x M.+ with values in M. d satisfying 
the following assumptions: 

• There exists 5q > such that 

(61) inf rnjti (|^|(x >7 ) + |^|(x,7)) l " !: ^i 00 +oo. 

• For every compact set K , 

(62) sup \h{(x,i) -^(x, 7)|^0. 

Then, for every sequence (x^)^^ ofR d , 
1 " 

— ^Vk sup ||$(/ i ?(x fc _ 1 , 7fc ))-$(^( 3 ; fc _ 1 ,7 fc ))||*(x fe _i) n ^±> oo 0. 

k=1 G [O,1] 

Proof. <1> has a compact support, therefore, we derive from (61) that 
there exists Mg > such that, for every |x| > M$ , 7 < <5o and 6 E [0, 1], 

*(/i?(x,7)) = 7)) = 0. 

Consider p 1— > w(p,§) = sup{r/ > 0, sup|.,,_ y | <r) |<3?(x) — <3?(y)| < p}. As <5 is 
uniformly continuous, w(p, <£) > for every p > 0. Thanks to (62), for every 
p > 0, there exists 5 P < 5o such that, for every 7 < <5 p , 9 £ [0, 1], 

sup \h l (x,-f)-h e 2 (x,-y)\<w(p,$). 
\x\<M So 

As 7^ fc 3L±f° Q ; there exists k p £ N such that 7^ < 5 p for k > k p . By using 
that i7 n n zL+f° +00, we deduce that 

1 n — 
limsup — sup ||$(h?(x fc _i, 7fe )) - ${h e 2 {xk-ulk))\\®{xk-i) < P^s<_ v 

n-»+oo -H n k=l 0e[O,l] 

where := sup{|^I / (x)|, |x| < M$ } < +00 since * is locally bounded. The 
result follows. □ 



PROOF of Proposition 4. We have to inspect successively the q € 
(1/2,1] and q£ [0,1/2] cases. 
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Case q £ (1/2, 1]. Let / € Cj^(R. d ). Decompose the infinitesimal generator 
as the sum of three terms defined by 

Aif(x) = (Vf,b)(x), A 2 f(x)=Tr(a*D 2 fa)(x), 



As/O) = J (f(x + K{x)y) - f(x) - (V/(x), K{x)y)l{\y\< l })-p:{dy). 

Set X kj i = X k __i + j k b(X k _i), X ki2 = X kA + ^y^a^k-^Uk and X k) 3 = 
X kt 2 + ^/lk^{X k -\)Z k . We then part the proof into three steps: 

Step 1. —E{f(X k>1 ) - f(X k _ 1 )/T k _ 1 } = A 1 f(X k _ 1 ) + R 1 ( lk ,X k „ 1 ) 
Ik 

1 n - n-oo 

with -Fr"l2' r ikRi(ik,x k _ 1 ) "-^o. 

Step 2. -E{/(X fei2 )-/(A ? fcil )|^_ 1 } = ^ 2 /(A > fc _ 1 ) + J R 2 ( 7fc ,A > fc _ 1 ) 
7fc 

1 ™ - n— oo 

With -7T^l^kR2{lk^X k _i) "-^ 0. 

Step 3. -E{f(X k )- /(4 2 )|fn} = A 3 /(A > fc _ 1 ) + J R 3 (7fc,^-i) 
7fc 



1 ™ - n— oo 

with — ^r? fc i? 3 (7fc,X fc _i) "-^0. 



The combination of the three steps yields Proposition 4. We refer to Proposi- 
tion 4 of [15] for steps 1 and 2 and focus on the last step where the specificity 
of our jump Levy setting appears. Since X k _i is J^-i-measurable and Z k , 
U k and T k -\ are independent, we have 

nf(Xk,2 + K(AVl)^)l^fc-l} = Q-yJ(Xk-l), 

where Qjf(x) = / E{f(S x ~, u + K(x)Z 1 )}P Ul (du), 

JR d 

with S x ^ tU = x + jb(x) + ^/7~ct(x)m. Set Vt = S x ^ >u + n(x)Z t . Applying Ito's 
formula to (f(V t )) t >o yields 



/G4) = /(S W )+ / <V/(V,-),K(x)dy s ) 

(63) _ J ° 

+ H f (S xa>u + K(x)Z s -,x,AZ a ) 

0<s<t 

where, H f (z,x,y) = f(z + K(x)y) - f(z) - (V f(z), K(x)y)l { \ y \< 1} . 

(64) 
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The process (/q (V/(Vg-), k(x) dY s )) is a true martingale since V/ is bounded. 
The compensation formula and a change of variable yield 

E{/(#r,7,u + k(x)Z 7 )} 

= /(^,7,u) +7 E |^ du y n(dy)H f (S x ^, u + k{x)Z v1 , x,y)j. 

Since A3 /(a;) = / ir(dy)H^ (x,x,y) = K{J^dv / ir(dy)H^ (x,x,y)}, it follows 
from the previous inequality that 

— E{/(X fc ) - = A 3 /(X fe -i) + i? 3 (7fc,^-i), 

7fe 

where, 

#3(7, a) = /e( / du / vr(dy)A^(5 w + K(x) 



with AH? (zi, Z2,x,y) = (zi,x,y) — H* (z2,x,y). We upper-bound i? 3 by 
two terms: -R31 and i? 3 2 that are associated to the small and big jumps 
components of (Z t ), namely, 

#3,1(7*3) 

= [ [ dv [ Tr(dy)lE\AH f (S Xt y tU + K(x)Z in ,x,x,y)\F Ul (du), 
J Jo J{\y\<i} 

#3,2(7^) 

= [ [ dv [ Tr(dy)K\AH f (S xn ^ u + K(x)Z in ,x,x,y)\F Ul (du). 
J Jo J{\y\>i} 

We study successively i?3 5 i and From Taylor's formula, we have for 

every y such that \y\ < 1 

\AH f (S X)J)U + K(x)Z vy ,x,x,y)\ < ±R(Z,-f,x,u,v,y)\K,(x)y\ 2 , 

where 

R(Z,j,x,u,v,y) 

= sup \\D 2 f(S xn!U + K(x)(Z vy + 6y))-D 2 f(x + 6K(x)y)\\. 
6»e[o,i] 

By setting $ = D 2 /, \P(x) = ||K(x)|| 2 |y| 2 , 

hi(x,^) = S x ^ tU + K(x)(Zi, 7 + 0y) and h^x , = x + 8 n(x)y , 

we want to show that the assumptions of Lemma 6 are a.s. fulfilled for every 
fixed u, v and y. 



(65) 
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First, since k(x) ' ' = + o(|x|), there exists a continuous function e such 

that k(x) = \x\e(x) and e{x) ' - — > 0. Therefore, as b and a have sublinear 
growth, one checks that there exist some positive real constants C\ and C 2 
such that 

|5x l7l « + k(x)(Z v1 + 9y)\ > \x\{l - 7C1 - {\Z V1 \ + \y\)\e{x)\) - C 2 , 
\x + 9k(x)u\ >\x\(1 — |e(x)||y|). 

Let So be a positive number such that 1 — 5qC\ > 0. Since {Zt) is locally 

bounded (as a cadlag process) and e(x) x - — >° 0, there exists a.s. M > 0, 

inf (l- 7 d- (\Z^\ + \y\)\s(x)\) >0. 

|x|>M,7e[0,5 ] 

It follows that a.s., 

inf (|fc?|(x, 7 ) + |^|(x,7)) +00. 
ee[o,i],7e[o,5o] 

Second, let if be a compact set of M. d . We check that (62) holds. We have 

sup \hx(x,j) -h e 2 {x)\ 
xeK,ee[o,i] 

(66) 

< sup(7|6(x)| + ^/7~||ct(x)|||'u| + ||k(x)|| \Z V1 \) — >0 a.s. 

x£K 

because b, a, k are locally bounded and lim^o Zt = a.s. Thus, by Lemma 
6, for any sequence (x k ) k& ^ of M. d , for every (u,v,y) Gl^x [0, 1] x B^O, 1), 

(67) — ^r] k AH f (S Xk _ ulktU + K(x k -i)Z vlk ,x k - 1 ,x k _ 1 ,y) ,1 -^?0 a.s. 
Hn k=i 

Now, since V/ and D 2 f are bounded, we derive from Taylor's formula that, 
for every z\, z 2 € W 1 , 

\frft \ ufl Mi . J 2 II V /IIoo||k(x)|||2/|1 { | |< 1} , 

Then, for every q € (1/2, 1], 

(68) |Aj^(S w + ^x)^ 7 ,x,:E,y)|l {M < 1} < C|Kx)|| 2 %| 2 n {M < 1} , 

where C = 2max(||V/||oo) ||-D 2 /||oo)- Therefore, by assumption (H 2 ), we fi- 
nally derive from (67), (68) and from the Lebesgue dominated convergence 
theorem that 

(69) 

1 

if sup — ^2r] k \\K,(x k -i)\\ 2q < 00. 

neN tin k=l 
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We apply this result to (x k ) = (Xk). By (58), sup ngN ^n(IM| 2l? ) < o° a.s. 
Hence, it follows that 1/H n J2k=i VkR3,i(lk, -Xfc-i) a.s. 
Now, let us focus on i?3 ; 2- Set A/ (z±,Z2) = f(z\) — f(zz). Then, 

1*3,2(1, x) = / E ( / dv [ Tr(dy)Af(x,S xn>u + K(x)Z vl )\¥ Ul (du) 
J Uo J{\y\>i} J 



+ [e( f 1 dv f ir(dy) 
J Uo Alyl>i> 



'{|i/l>i} 

x Af(S XjljU + k(x)(Z v1 + y),x + K,(x)y)^¥ Ul (du). 

One proceeds as before. By using Lemma 6, one begins by showing that, for 
any sequence (xfc)fceN , for every (u,v,y) 6 [0, 1] x M. d x Bd(0, l) c , a.s., 

1 ™ n^oo 
7T X! 1^/(^-1. S in,7b« + K fe-l) 2 «7fc ) ™ alld 

1 

— ^2r) k Af(S Xk _ ulk , u + «:(xfc-i)(^ 7fe +y),x k -i + K(x k -i)y) 0. 
Hn k=i 

By the dominated convergence theorem [which can be applied because 
Tt(\y\ > 1) < 00 and / is bounded], we deduce that, for any sequence (xk)i< 

1 n 



— VkR3,2(lk,Xk-l) "-^ a.s. 

-tin. , 



fc=l 

This completes the proof of Step 3 when q£ (1,2]. 

Case q < 1/2. The reader can check that the assumption q £ (1/2, 1] is 
used only once: when we want to apply the dominated convergence theorem 
for i?3 5 i [see (68)]. Since inequality (68) is not true when q < 1/2, we need 
to decompose the infinitesimal generator in a slightly different way: 

A 1 f(x) = (Vf,b°)(x), 
A 2 f(x) = Tiia*D 2 fa)(x), 

Mf(x) = [ (f(x + K(x)y) - f(x))v(dy). 



Note that this decomposition is only possible when q < 1/2. That means 
that with the notation (40), we decompose AXj. with h = and inspect 
the three induced steps. We do not go into further details since the proof is 
similar to the case q > 1/2. □ 
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Proposition 5. Assume that the assumptions of Proposition 3 are ful- 
filled. Then, 

lim — £ ^E{f(X k ) - /(Xfc-OpVl} = a.s. 

Proof. We do not detail the proof of this proposition which is an adap- 
tation of Proposition 3 in [15]. □ 

5. Proof of the main theorems for Schemes (B) and (C). The aim of 

this section is to give a general idea of the proof for Schemes (B) and (C) 
and to overcome the main difficulties induced by the approximation of the 
jump component. For Scheme (A), main theorems have been proven in two 
successive steps. First, we focused on tightness results (Proposition 1) and 
then proved that every weak limiting distribution is invariant for (Xt)t>o 
(Proposition 3). We follow the same process for Schemes (B) and (C). We 
will successively explain for both schemes why Proposition 1 and Proposition 
3 remain valid. 

5.1. Almost sure tightness of u^(u,dx) and u^(u,dx). The tightness 
result for Schemes (B) and (C) is strictly identical to Proposition 1 [in 
particular, assumption (9) is not necessary for tightness]. Looking carefully 
into the proof of this theorem for Scheme (A) shows that the properties of 
the jumps that we use are the following: the control of the moments of the 
jump components (Lemma 2) which is fundamental for Proposition 2, and 
independence between (Y n ) n& ^, (N n ) n ^n and (U n ) ne ^. We show in Lemma 
7 below that the controls of Lemma 2 hold true for the moments of the jump 
components of Schemes (B) and (C). Then, since Scheme (B) satisfies the 
independence properties, Proposition 1 follows in this case. In Scheme (C), 
(Y^p) n6 N and (Nn)n£N are no longer independent. It raises several technical 
difficulties in the proof of Proposition 2 in case p < 1, but the process of 
the proof is the same. So, we only state a variant of Lemma 2 (see [21] for 
details). 

Lemma 7. Let Tq be a positive number and T n = inf{s > 0, |AZ™| > 0}. 

(i) Let p>0 such that (Hp) holds. Then, for every t < To and h > 0, 

n\Nt T A 2p }<t I \y\ 2p v(dy) i f P >0. 

J\y\>h 

(ii) Let t be an (JFt)- stopping time and q 6 [0, 1] such that (Hq) holds. Set 
D n = {y,\y\ G (u n ,h]} and Y t h ' n = T,o<s<t AY s hl {AYheDh} ~ tf D ny^{dy). 
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Then, 

2q 

E" ' 



Y t ^ + (tAr)[ y7r(dy) q \<t[ \y\ 2(! n{dy), i/ge [0,1/2], 
E{ |yM|2« } < C i J \ y \*H{dy), ifqe (1/2,1]. 

(iii) Letp> 1 such that (K^) holds. Set Zf = Zf — t /{u>u V^idy). Then, 
there exists r\ > 1 sitc/i i/iaf, /or ewer?/ T > 0, /or ener?/ e > 0, i/iere exisis 
Ce,T ,p > 0, no £ N s-uc/t i/iai, /or ewery t > To and n > uq, 

E{|Zlf p }<*(/ M 2p 7r(^) + e) +C £!T0)P i" 

and 

E{|ir A T"| 2p }<i(/ |y| 2 M^) + e) 
Proof. The proof is left to the reader. □ 

Remark 9. In (iii), the control is only valid for n sufficiently large but 
that does not make any problem since (R a . P ) just needs to be valid for 
sufficiently large n. 



5.2. Identification of the limit of (u B ) n ^ and (u^) n& ^. The theorem 
which is obtained for (u B ) n ^ and (z^jneN is strictly identical to Proposition 
3 under the additional condition (9) for Scheme (C). We recall that the proof 
of Proposition 3 is based on two steps: Propositions 4 and 5. Proposition 5 is 
still valid without additional difficulties. However, the proof of the analogous 
result to Proposition 4 raises some new difficulties. Denote by A k,B and A k,c 
the operators on Cf (R d ) with values in Cb(R d ,IR) defined by 

A k > B f(x) = (Vf,b)(x) + \Ti{a*D 2 fa)(x) + f & Xx,x,y)*(dy), 

A k ' c f(x) = A k > B f(x) - (1 - a k ( lk )) ( Hf(x, x, y)vr(dy), 

J{\y\>u k } 

where a k (t) = ^f^^f ■ " : A f - A k > B /" and 11 A f - A k > c /" can be viewed 
as the principal part of the weak error induced by the approximation in 
Schemes (B) and (C) [A k,B is the infinitesimal generator of (X k ), where 
(X k ) is solution to the SDE (1) driven by (Z k ) instead of (Z t )]. Thus, one 
may expect that this error be negligible in the sense of our problem. This is 
the aim of Lemma 8. 
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Lemma 8. Assume (Hq). Let (x k ) k ^ be a sequence such that 

1 n 

(71) sup — Y,Vk\\< x k-i)\\ 2q <oo. 



n>l H n 



Then, for every function f £ (M d ,R), 



and i/ 7r(D n )j n "zli? q, 



1 n 

lim — J2 m(Af(x k ^) - A k ' B f(x k ^)) = 



n 

lim — J2 »7fc(A/(x fc _i) - A*- c /(x fc _i)) = 0. 



fe=l 



PROOF. Note that A k ' B f(x) -Af(x) = fn y \ <Uk \ (x, x, y)ir{dy). When 
q > 1/2, we deduce from Taylor's formula and the boundedness of V/ and 
D 2 f that there exists C q > such that 

|^ / (a; ) x,y)|l { | y |< ttfc} <C F ,||K(x)|| a «|y| 2 «l { | 1/ |< ttik} . 
When q < 1/2, since / is a 2(/-H61der function, 

|^(x,x,y)|i {H <„ fc} < [fUMx)\\ 2q \y\ 2 n {muk} 

+ SUp |V/(x)| • ||«(X)|| |j/|l{|y|<u fc }- 
xSsupp / 

By setting v Kq = J { \ y \ <Uk} \y\ 2q ir{dy), we have 



Since v k>a k -^f for every a > q under assumption (Hq), the first result 
follows from (71). One deduces the second inequality by checking that 

\A k > B f(x) - A k > c f(x)\ < Cir{D k ) lk {l + ||k(x)|| 2 "). □ 



Set 



T>B,k Ui yB \ E if( X k ) ~ /(^fc-l) 1-^-1 i \h,B f( v~B \ 

Ik 

t? c > k (~ y-c \ E if(Xk) - /(^"fcLi)|-^i} A k,Cf(yC \ 
K 3 Uk, x k -i) = A J{ x k-i)- 

Ik 

The rest of the proof then amounts to proving that 



1 n 
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and 



1 



lim TT £ VkRrilk, Xg. x ) = 0. 



- H n fe=1 

We do not detail this proof based on the same approach as the proof of 
Proposition 4 (see [21] for more details). However, we want to derive the 
main difficulties from the proof. For Scheme (B), one deduces from the Ito 
formula that 

\R k 3 ' B (7,x)\ < J J\v J 7r(dy)E\AHf(S X)1)U + K(x)Z^,x,x,y)\P Ul (du). 

The right-hand term can be written i?^ fc (7, x) + R^ 2 k (l, x), where R^'i(^y, x) 
[resp. Rz2{li x )\ 1S simply derived from ^3^(7, x) [resp. ^3,2(7, x)], defined 
in the proof of Proposition 3, by replacing Z with Z k . We focus on R^i ■ One 
observes that the controls (65) and (66) used for R31 no longer work since 
the jump component depends on n. An idea is to use the Skorokhod repre- 
sentation theorem (see, e.g., [24]) in order to replace (Z k ) by a uniformly 
controllable sequence. 

Lemma 9. There exist a sequence of cdlg processes (Z n ) and a cdlg 

process Z such that Z n = Z n for every n > 1, Z = Z and Z n — > Z a.s. for 
the Skorokhod topology. In particular, 

sup sup I Z™ I < +00 VT > and 

neN0<s<T 

(72) 

limsup sup |Z™|=0 a.s. 

+0 0<s<7 

Proof. Z n converges locally uniformly in L? toward Z, hence, in distri- 
bution for the Skorokhod (Polish) topology. Thanks to the Skorokhod rep- 

~ ~ ~ c ~ c 

resentation theorem, there exists (Z n ) ng pj and Z with Z n = Z n and Z = Z 

such that Z n tends a.s. toward Z for Skorokhod topology. The assertion 

(72) easily follows from the continuity of <y 1 — > || ol || gup and ol 1 — > o;(0) for the 

Skorokhod topology. □ 

Since R%'i only depends on the law of Z n , one can replace Z n with 

Z n . Then, we use (72) as an alternative to the local boundedness and the 
continuity at t = of (Z t ) needed in (65) and (66) respectively. A result 
analogous to (67) follows. The idea is the same for R32 ■ 

Finally, for Scheme (C), the result essentially follows from the following 
remark: 

sup \Z™ AT „ \ < sup \Z™\. 

0<s<t 0<s<t 
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This means that the remainders in Scheme (C) are easier to control than 
those of Scheme (B). For more details, we refer to [21]. 

6. A theoretical application. The "classical" a.s. CLT due to Brosamler 
[6] and Schatte [26] is the following result. Let (U n ) n ^ be a sequence of 
i.i.d. random variables with values in R rf such that KUi = and T i jj 1 = Id- 
Then, 

1 n 1 

P-a.s V-<5 i//r _ fe =>jV(0,I d ). 

Inn £~t k VvfeEfci^ 

k=l 

This result is obviously connected with the central limit theorem which 
expresses the fact that every square-integrable centered random variable 
is in the domain of normal attraction of the normal law. When the square- 
integrability no longer holds, Berkes, Horvath and Khoshnevisan [3] obtained 
an extension of this result connected with the nonsquare-integrable attrac- 
tive laws which are stable laws [with index a € (0, 2)]. We are going to show 
that we can deduce this extension from Theorem 2. 

Let c be a positive number and denote by (Z"' c )t>o a symmetrical one- 
dimensional a-stable process such that the characteristic function (j) of Z^ ,c 
satisfies 4>{u) = e _p ' u ' a , where p = 2c / +O ° U~ a sinydy. Consider a sequence 
(Vn)neN of symmetrical i.i.d. random variables such that, for x > 0, 

P(y x > x ) = — + 5{x){x^ a {lnx)^) 

(73) 

with 7 > and S(x) — ► 0. 
By a result of Gnedenko and Kolmogorov (see [10]), we know that 

VjJ Wn ^ 7 a, c 

Then, the following a.s. CLT holds: 

Theorem 4. Let {j]k)k£M be a nonincreasing sequence with infinite sum 
such that (A:r/fe)fc g N is nonincreasing and set v = C(Z^ ,C ). Then, if r y>—, 
a.s., 

i A m 

jr 2^ ^°(v 1 +-+v k )/ky^ =^ v - 
Un k=i 



In particular, 



RECURSIVE COMPUTATION 



41 



In order to prove this theorem, we first need an almost sure invariance 
principle due to Stout (see [27] or [3]). 

Proposition 6. Let (V n ) n >i and (Cn)n>i be sequences of i.i.d. random 
variables such that d = Z^' c and V\ is defined as above. Then, if^> — , there 
exists a probability space (fl,^F,¥) and sequences of i.i.d. random variables 
{Vn)n>i and (Cn)n>i such that V\ = Vi, Ci = Ci and 

n n / 1 \ 

( 74 ) E & - E V* n ^-°° oin^ilnnm a.s. Vp G 0, 7 - - . 
i=l i=l \ CtJ 

Proof of Theorem 4. First, we assume that Vi = Ci = Z?' c . Set 

„ Cl H 1" Cn+l w . „ 

S n = — ; -—, Vn > 0. 

(n+l)V« 

i_ 

One easily checks that S n+1 = S n - ±7 n+ iS n + ^ +1 ( n+2 + Rn+i with -y n = 
and -R n +i = 0(j 2 +1 \S n \). The idea of the proof is to compare (S n ) n >o 
with the exact Euler scheme with initial value £j and step sequence ( 7n ) 
associated with the SDE (E Q]C ) defined by dXt = ——X t - dt + dZ^' c . Since 
(Z"' c ) t >o is a self-similar process with index 1/a (see, e.g., [25]), its Euler 
scheme can be written 

1 _ I 

X = Ci and X n+i = X n j n +iX n + 7^ +1 Cn+2- 

a 

As an Ornstein-Uhlenbeck process driven by a symmetric stable law, (Xt) 
admits a unique invariant measure v and v = C(Z"' C ) (see [25], page 188). 
Since k is bounded, assumptions of Theorem 2 are clearly fulfilled with 
V{x) = 1 + x 2 , a = l and for any p € (0, a/2) and q £ (a/2, 1). (In the rest of 
the paper the initial value of the Euler scheme is supposed to be constant, 
but it is obvious that Theorem 2 is still true when X$ is a random variable 
satisfying E{|A > o| 2p } < +oo.) Hence, it follows from Theorem 2 that 



a.s. 

71 k = l 



Then, by using that \f(Sk) — f{Xk)\ < C(\Sk — X^\ A 1) for every Lipschitz 
bounded function /, one easily checks that Theorem 4 holds with V\ = Ci if 

(76) A n :=S n -X n n ^°0 a.s. 

Let us show (76). One first checks that 

A = and A n = f 1- 1 )A n _i + #» Vn > 1. 
V a(n + l)J 
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Setting fco = inf{fc > 0, k — 1/a > 0}, we deduce that, for every n > kg + 1, 

Al n n f 

A n = -^ + - ^ c kRk withc n = 1 



One observes that 



cxp 



jLK l -=)M=jLM*))^^- 



Then, Afe /c n — > a.s. Hence, (76) holds if we check that l/c n J2k=k +i c kRk — 
a.s. First, if a > 1, as £i is integrable and Rk = 0(Sk-i/(k + l) 2 ), we have 



fe>l k> 

We deduce that J2k>i \Rk\ < +oo a.s. Since c n n zL±f° +00, we derive from 
Kronecker's lemma that 

i n 

— 2^ Cfc/tfc — > U A n — > U a.s. it a > 1. 

Cn fc=fc +i 

Second, if a < 1, Ci has a moment of order for every < a. It follows from 
inequality (25) that 

n\Rkf} < C - ; AzT^7 ^{ICl| e } < 



(ifc + 1)20+0/^1^1 / - (fc + 1 )0(2+l/a)-l " 

Therefore, if 6 satisfies 0(2 + ±)-l>l, that is, if ^ < # < a, we have 
X)a,>i \Rkf < +oo a.s. Hence, by inequality (25) and Kronecker's lemma, it 
follows that 



1 - 

— ^2 Ck ^ k 
0x1 k=k +i 



1 n 

4 E 4\R k f n ^0 a.s. 
n fc=fc +i 



and the theorem is proved when V\=Q\. Now, consider a sequence (V n ) n >o 
of i.i.d. symmetric random variables satisfying (73). Since Theorem 4 is 
true for (Cn)n>i, it is also true for every sequence (Cn)n>i of i.i.d. random 

variables satisfying Ci = Ci- By taking (Cn)n>i such that Proposition 6 holds, 
we derive from (74) that there exists a sequence of i.i.d. random variables 
{Vn)n>i such that V\ = V\ and 
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( 77 ) -ft m$ i vi+...+v k )/ki/« n ^°° v a - s - 

nn k=l 

As (V n )n>i and (V n ) n >i are sequences of i.i.d. random variables such that 
Vi = Vi, (77) is also true for (V n ) n >i. □ 

7. Simulations. 



Example 4. Denote by (Zt)t>o a Cauchy process with parameter 1 
[with Levy measure defined by n(dy) = 1/y 2 dy] and consider the Ornstein- 
Uhlenbeck process solution to dX t = —X t ~ dt + dZ t corresponding to (Ei t \) 
defined in the previous subsection. The unique invariant measure of (Xt)t>o 
is the Cauchy law (see [25], page 188) and the assumptions of Theorem 2 
are fulfilled with V(x) = 1 + x 2 , a = 1 and every p £ (0, 1/2) and q £ (1/2, 1). 
Therefore, 

for every / satisfying / = 0(|x| 1//2_e ) with e > 0. In Figures 1, 2 and 3, 




-20 -15 -10 -5 5 10 15 



Fig. 1. Scheme (A), t = 12.5. 
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— Theor, density 



-20 




Fig. 3. Scheme (C), t = 16.4. 
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Fig. 4. Scheme (A). 



one compares the theoretical density of the invariant measure with the den- 
sity obtained by convolution of each of the empirical measures by a Gaus- 
sian kernel for N = 5.10 4 . We choose % = 7 n = 1/y/n, u n = [so that 
7r(D n )-y n — > 0] and t indicates the CPU time. In order to have a more precise 
idea of the differences between the three Euler schemes, we simulate and 
represent on Figures 4, 5 and 6 the sequence {v n {f)) with f(x) = \ x\ 0A , for 
several choices of polynomial steps. We set j n = r\ n = 1/n 9 and u n = 7 n (resp. 
u n = y^y ) for Scheme (B) [resp. for Scheme (C)]. We observe that, among 
the tested steps, the best rate seems to be obtained for 6 = 0.3. Notably, in 
Schemes (B) and (C), we see that, on the one hand, if the step decreases 
too slowly (e.g., = 0.7), so is the stabilization and, on the other hand, 
if the steps decreases too fast (e.g., when 9 = 0.1), there are not sufficient 
variations to correct the error. 

Remark 10. In [20] we study the rate of convergence of these procedures 
in terms of steps, weights and truncation thresholds. This enlightens these 
first numerical illustrations. 

Example 5. Now we deal with the following SDE: 
dX t = (l-X t -)dt-X t -dZ t , 
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Fig. 6. Scheme (C). 
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2.5 



1.5 



0.5 



Scheme (B) 

Scheme (C) 



0.1 0,2 0-3 0.4 0-5 0.6 0.7 0,8 
Fig. 7. Approximated density, N — 10 e '. 



0.9 



where (Zt)t>o is a drift-free subordinator with Levy measure tt defined by 

x /3/2,l/2(y) , 

K{dy) = 2 dy, 

where f a ^ is the density function of the f3(a, ^-distribution. This SDE mod- 
els the dust generated by a particular EFC process (see Introduction) whose 
sudden dislocations do not create dust, having parameters (according to the 
notation of [2]): 



Ck = 0, 



1. 



Vcoa 9 {dy) = h/2,i/2(y)dy. 



One checks that (S 11:L /2) is satisfied with V(x) = 1 + x 2 . However, we do 
not have k(x) =o(\x\), but since supp(7r) is restrained to [0,1] without sin- 
gularities in and 1, we are able to show that assumption k(x) = o{x) is no 
longer necessary in this case. In Figure 7 we represent the approximation of 
the invariant measure obtained for Schemes (B) and (C) [we are not able to 
simulate Scheme (A) in that case]. 



Acknowledgments. Thanks to Gilles Pages for extensive discussions and 
suggestions. 
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